A  FOS  R -  1   13  7 


NEW  YORK  UNIVERSITY 

INSmrUTE  OF  MATHEMATICAL  SCKNCIS 

LIBRARY 

(«  WaAinglon  PUce,  New  Yot4(  3,  M  Y, 

OCT  23  10*^1 


^^E  ET  PR^^  NEW    YORK    UNIVERSITY 

-    KrT\/5 

w   1    N   11      ^  Institute  of  Mathematical  Sciences 

'V/s  _  _  _  ^  V  4^  Division  of  Electromagnetic  Research 


RESEARCH    REPORT    No.    HT-9 


Stochastic  Models  for  Many-Body  Systems 
I.  Infinite  Systems  in  Thermal  Equilibrium 


ROBERT  H.  KRAICHNAN 


Contract    No.    AF    49(638)-341 
Project    No.    47501 
JULY,     1961 


NEW  YORK  UNIVERSITY 

Institute  of  Mathematical  Sciences 

Division  of  Electromagnetic  Research 

Research  Report  No.  HT-9 

July,  1961 


STOCHASTIC 

MODELS  FOR 

MANY -BODY  SYSTEMS 

I.  INFINITE 

SYSTEMS 

IN 

THERMAL 

EQUILIBRIUM 

Robert  H 

Kraichnan 

li  ^^^.U. 


Robert   H.    Kraichnan 


^TH^     "^^ 


Morris  Kline;,  Director 


Qualified  requestors  may  obtain  copies  of  this  report  from  the 
Arlington  Hall  Station,  Arlington  12,  Virginia.   Department  of 
Defense  contractors  must  be  established  for  ASTIA  services,  or 
have  their  'need  to  know'  certified  by  the  cognizant  military 
agency  of  their  project.   Reproduction  in  whole  or  in  part  is 
permitted  for  any  purpose  of  the  U.S.  Government. 


ABSTRACT 

Some  model  Hamiltonians  are  proposed  for  quantum-mechanical  many- 
body  systems  with  pair  forces.   In  the  case  of  an  infinite  system  in 
thermal  equilibrium,,  they  lead  to  temperature-domain  propagator  expan- 
sions which  are  expressible  by  closed,  formally  exact  equations.   The  ex- 
pansions are  identical  with  infinite  subclasses  of  terms  from  the  propagator 
expansion  for  the  true  many-body  problem.   The  two  principal  models 
introduced  correspond,  respectively,  to  ring  and  ladder  summations  from 
the  true  propagator  expansion,,  but  augmented  by  infinite  classes  of  self- 
energy  corrections.   The  latter  are  expected  to  yield  damping  of  single- 
particle  excitations.   The  eigenvalues  of  the  ring  and  ladder  model 
Hamiltonians  are  real,  and  they  are  boionded  from  below  if  the  pair 
potential  obeys  certain  conditions.   The  models  are  formulated  for  fer- 
mions,  bosons,  and  distinguishable  particles.   In  addition  to  the  ring 
and  ladder  models,  two  simpler  types  are  discussed,  one  of  which  yields 
the  Hartree-Fock  approximation  to  the  true  problem.   A  novel  feature  of 
all  the  model  Hamiltonians  (except  the  Hartree-Fock)  is  that  they  contain 
an  infinite  number  of  parameters  whose  phases  are  fixed  by  random  choices. 
Explicit  closed  expressions  are  obtained  for  the  Helmholtz  free  energy 
of  all  the  models  in  the  classical  limit. 
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1 .    INTRODUCTION 

A  difficulty  in  calcxilating  the  statistical  behavior  of  many-body 
systems  is  that  some  quantities  of  interest  may  not  have  convergent  per- 
turbation expansions.   If  the  system  is  infinite  in  size,  this  can  occur 
even  when  the  density  of  particles  is  very  low  and  the  interparticle 
forces  are  very  weak.   The  situation  is  already  present  in  the  classical 
theory.   Consider  a  gas  of  classical  particles  which  interact  by  a  re- 
pulsive, short-range  pair-potential  XV(x).   If  the  pressure,  expressed 
as  a  function  of  density  and  temperature,  is  expanded  in  powers  of  the 
strength  parameter  X,  the  expansion  has  zero  radius  of  convergence. 
This  suggests  that  a  similar  situation  may  exist  in  the  quantum-mechanical 
case. 

Non- convergence  of  interaction-strength  expansions  is  not  necessarily 
a  disaster.   In  the  classical  example  just  cited,  the  expansion  almost 
certainly  is  asymptotic  about  X  =  0,    and  we  may  hope  that  this  is  true 
also  in  some  quantum-mechanical  cases.   However,  many  physical  problems 
of  interest  do  not  exhibit  weak  interactions.  Moreover,  certain  proper- 
ties of  a  quantum-mechanical  many-body  system  may  not  have  even  asymptotic 
expansions  as  power-series  in  X.     An  example  is  the  one-particle  momentum 
distribution  (?(k)   normalization:  /  cr(k)d  k  =  1  .   For  an  infinite  system, 
a  finite  change  in  (?(k)  from  its  form  for  lincoupled  particles  means  that 
an  infinite  number  of  particles  are  displaced  from  the  momentum  levels 
they  wo^old  occupy  if  they  were  not  coupled.   In  order  to  form  such  a  state 
from  the  uncoupled  state,  the  interaction  Hamiltonian  must  act  an  infinite 
number  of  times;  that  is,  infinite  orders  of  perturbation  theory  are  in- 
volved. We  expect  on  physical  grounds  that  the  change  in  (^(k)  goes  to 


zero  as  X  does.   However,  we  cannot  presume  that  it  must  go  to  zero  as 
some  integral  power  of  X. 

In  recent  years,  several  formalisms  for  handling  perturbation  expan- 
sions in  the  quantum-mechanical  many-body  problem  have  been  proposed  which 
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are  related  to  methods  previously  used  in  quantum  electrodynainics . 

They  produce  great  simplifications  in  manipulations  and  permit  one  to 
carry  out  various  formal  summations  of  infinite  classes  of  terms  from 
the  perturbation  expansions.   These  formalisms  are  a  natural  choice  for 
investigating  quantities,  such  as  (3'(k),  which  may  not  have  convergent  or 
asymptotic  expansions.   However,  it  can  be  difficult  to  know  in  advance 
what  summations  should  be  carried  out  in  given  cases.   With  sufficient 
ingenuity,  it  is  possible  to  sum  certain  divergent  series  so  as  to  obtain 
almost  any  answer  whatever.   It  is  difficult  to  guess  in  advance  whether 
adding  further  infinite  classes  of  terms  to  a-  known  expansion  will  improve 
the  answer  or  make  it  worse . 

In  the  present  paper  we  shall  describe  a  procedure  intended  to  pick 
out  perturbation-term  summations  for  which  certain  characteristics  are 
predictable  in  advance.   We  shall  formulate  model  Hamiltonians  such  that 
the  complete  perturbation  expansions  to  which  they  lead  are  formally  iden- 
tical with  certain  infinite  subclasses  of  terms  from  the  corresponding 
perturbation  expansion  for  the  true  Hamiltonian.   The  model  Hamiltonians 
are  Hermltian,  and  have  eigenvalues  which  are  bounded  from  below  if  the 
pair  potential  obeys  certain  restrictions.   The  predictable  characteris- 
tics of  the  model  solutions  are  those  which  follow  from  these  properties. 


In  common  with  the  true  many-body  Hamiltoniaxi,  our  models  (vrLth  one 
exception)  are  not  dlagonalizable  by  known  means.   However,  they  are  sol- 
uble in  the  sense  that  they  yield  formally  closed  integral  equations  for 
the  propagators  that  determine  the  mean  energies,  meaji  occupation  numbers, 
etc.  which  are  of  statistical-mechanical  interest.  A  novel  feature  of 
the  models  is  that  they  contain  infinite  niimbers  of  parameters  whose 
values  are  chosen  at  random.   We  shall  therefore  call  them  stochastic 
models.   The  random  parameters  will  be  described  in  Sees.  2  and  5- 

We  shall  present  two  principal  types  of  stochastic  models,  ladder 
and  ring.   They  correspond,  respectively,  to  summations  of  familiar  in- 
finite classes  of  ladder  or  ring  diagrams  from  the  perturbation  series 
for  the  true  Hamiltonian.  At  the  same  time,  however,  they  include  cer- 
tain infinite  classes  of  self-energy  corrections  to  these  diagrams.   The 
corrections  are  of  a  type  expected  to  contribute  to  the  damping  of  ele- 
mentary excitations.   In  addition  to  the  ladder  and  ring  models,  we  shall 
introduce  two  simpler  types.   One  yields  the  Hartree-Fock  approximation 
to  the  true  problem  (and  contains  no  random  parameters).   The  other  also 
includes  the  Hartree-Fock  diagrams,  but  with  iterated  self-energy  correc- 
tions.  The  eigenvalues  of  this  last  model  are  not  bounded  from  below 
for  any  pair  potential,  and  its  validity  therefore  is  quite  doubtful. 

The  models  which  are  described  in  the  present  paper  yield  closed 
equations  only  for  infinite  systems.   We  shall  develop  each  type  of  model 
in  two  forms:  for  indistinguishable  particles  (fermions  and  bosons)  and 
for  distinguishable  particles.   The  analytical  treatment  will  begin  with 
the  distinguishable  particle  models.   They  admit  more  immediate  physical 
interpretations.   We  shall  apply  to  them  the  Ursell-Mayer  irreducible 


-  k 


cluster  expansion  method  and  thereby  obtain  an  explicit  closed  expression 
for  the  Helmholtz  free  energy  of  each  model  in  the  classical  limit.   For 

the  fermion  and  boson  models,  we  shall  use  a  temperature-domain  propaga- 

2 
tor  formalism  of  the  type  originated  by  Matsubara  and  developed  further 

by  Fradkin,    Abrikosov  et  al.,    Luttinger  and  Ward,    and  others.   The 
distinguishable  and  indistinguishable  particle  models  turn  out  to  give 
formally  identical  thermodynamics  in  the  classical  limit.   Thus,  our  classi- 
cal results  for  the  Helmlioltz  free  energy  provide  some  insights  into  the 
behavior  of  the  fermion  and  boson  models. 

In  the  paper  which  follows,  we  develop  more  general  models  which 
yield  closed  equations  whatever  the  size  of  the  system.   We  shall  apply 
them  to  non-equilibrium  as  well  as  equilibrium  statistical  mechanics. 
For  an  infinite  system  in  equilibrium,  the  generalized  models  yield  the 
same  final  equations  as  the  models  of  the  present  paper.   However,  they 
provide  a  neater  and  more  satisfactory  derivation  of  these  equations. 
We  do  not  start  with  the  general  treatment  in  the  present  paper  because 
it  requires  a  more  elaborate  formalism  and  therefore  does  not  provide  as 
direct  an  introduction  to  the  use  of  stochastic  models. 

The  derivation  of  our  closed  model  equations  involves  a  deep-lying 
convergence  question  which  is  described  in  Sec.  5.1.   We  make  no  attempt 
to  answer  this  question  in  the  present  paper.   In  the  following  paper,  we 
offer  what  we  hope  is  a  satisfying,  although  non-rigorous,  resolution. 


2,   MODELS  FOR  DISTINGUISHABLE  PARTICLES 
2.1  Nature  of  the  Models 

Let  us  consider  a  system  of  N  similar  but  distinguishable  particles 
(of  unit  mass)  which  interact  through  a  pair  potential  V(x).  The  total 
Hamiltonian  may  be  written 


k4 

n 


H!  Pn^  +  i  XI!    V(^n  -  ^m^     (n,m  =  1,2,  .  .  .  ,N),    (2.1) 


where  /  ,   means  that  n  =  m  is  omibted  in  the  sum.  Here  p  and  x 
^ — '  -^n      n 

are  the  momentiom  and  position  of  the  nth  particle.   We  shall  adopt  the 
artifice  of  confining  the  system  in  a  cubical  cyclic  box  of  volume  fi.. 
That  is,  we  restrict  V(x)  to  the  form 

V(x)  =  ^^  V  exp(ik-x),  (2.2) 

k 

where  k  takes  all  the  values  allowed  by  cyclic  boundary  conditions  on 
the  walls  of  the  box,  and  we  require  that  the  Schrodinger  wave  fimction 
be  a  cyclic  function  of  the  coordinates  of  each  particle.   In  the  classi- 
cal case,  we  assiome  that  a  particle  which  exits  through  any  wall  of  the 
box  simxiltaneously  re-enters,  with  the  same  momentum,  through  the  oppo- 
site wall.   We  shall  eventually  be  interested  in  the  limit  N  — >  <»  , 
fi  -^  CO,  with  W/n  finite. 

We  require  that  V(x)  be  real  and  have  reflectional  symmetry.   These 
conditions  imply 


V(x)  =  V(-x),      Y^  =  V_^,  \=\    ■  (2.5) 

Except  where  we  specify  otherwise,  we  shall  assume  that  V(x)  is  a  smooth, 
bounded  function  such  that 

|Vj^|  =  0(n'^),      n^^    (all  k), 


|V^|  §  0(k"^),     k  ^«=  . 

In  particular,  this  implies  that     /  V(x)d  x   and    /  [v(x)]  d  x 
exist  for  Q   — >  <»,  where  the  integration  is  over  the  whole  box. 

We   shall   call    (2.1)    the   true   Hamiltonian,    and  refer  to  the   statis- 
tical-mechanical problem  associated  with  it   as   the   true  problem. 

Now  let  us  consider  model  Hamiltonians  of  the  form 

n  n,ra 

where  V  '  (x)  is  a  pair  potential  which  may  be  different  for  each  pair 
of  particles  n  and  m.  We  require  that  V  '  (x)  be  real,  and  we  replace 
(2.5)  by 

V"^"'(x)  =  v'"'"(-x).  (2.5) 

Let  us  write 

V"'"(x)  =  C  O"  exp(ik.x)  (2.6) 

k 

and  define  the  parameters  i£  ,  by 

^  ^n,m;k 


v,''^"  =  V,  <j^        ,  (2.7) 

Then  the  reality  condition  and  (2.5)  imply 

'^n^mjk  =  '^n,m;-k  '     ^n,m;k"  '^m,n;-k-  (2.8) 

The  true  problem,  of  course,  corresponds  to  jz^   -k  ~  "'"  ^'^^  ^'^^  ^'    ^' 
and  k.   In  the  models  with  which  we  shall  be  concerned  here,  all  the 
j^    ,   will  be  assigned  unit  modulus  (except  in  the  Hartree-Fock  model 
where  most  of  them  will  vanish) .   However,  for  each  triad  n,  m,  k  the 
phase  of  jZ^       will  be  assigned  by  a  random  choice,  subject  to  (2.8)  and 
to  additional  constraints  which  differ  for  each  model.   As  we  shall  see, 
the  models  so  produced  have  certain  properties  in  common  with  the  true 
Hamlltonian  but  lead  to  a  statistical  mechanics  that  can  be  expressed  in 
closed  form  in  the  limit  17  — >  <»  • 


2.2  Ladder  Model 

Let  us  specialize  <k  ,  to  the  form 

•^         '^n,m;k 

jZ^    ,  =  exp(-ik.d   ),       d    =  -d   ,  (2.9) 

n,m;k      ^     n,m  n,m     m,n 

the 
where^d    are  constant,  real  vectors.   This  clearly  satisfies  (2.8). 

Now,  for  each  pair  n,m  let  us  give  the  three  vector- components  of  d 
values  chosen  at  random  within  the  interval  (0,L)  where  L  =  fi   .   The 
choices  are  to  be  completely  independent  for  pairs  which  are  not  identi- 
cal.  In  the  X  representation  ,  we  have 


v'^'^Cx)  =  V(x  -  d^^^),  (2.10) 

which  permits  a  very  simple  interpretation  of  this  model:   The  pair  poten- 
tial has  the  same  shape  as  in  the  true  problem,  but  the  particles  now 
collide  with  ghosts  of  each  other,  displaced  by  the  randomly  chosen  vectors 
d    .   For  reasons  which  will  appear  later,  we  shall  call  this  the  ladder 
model. 

An  important  feature  of  the  ladder  model  is  immediately  apparent  from 
(2.10).   If  V(x)  is  non-negative  for  all  x,  then  V  '    (x)  also  has  this  pro- 
perty.  It  follows  that  in  this  case  the  expectation  of  H  in  any  quantiom- 
mechanical  state  is  non-negative,  or,  in  other  words,  that  all  the  eigen- 
values of  H  are  non-negative. 

At  this  point,  we  want  to  malce  as  clear  as  possible  the  precise  sense 

in  which  our  model  is  stochastic.   The  values  of  the  d  -    are  chosen  at 

n,m 

random  [subject  to  (2.9)] •   Once  chosen,  however,  they  are  fixed,  and  we 

work  thereafter  with  the  definite  Hamiltonian  embodying  these  values.   In 

particular,  the  same  choice  of  the  d    will  be  employed  for  every  member 

of  the  canonical  or  grand  canonical  ensemble  which  we  use  in  describing 

the  statistical  mechanics  of  the  system.   The  principal  deductions  we  shall 

make  about  the  statistical  mechanics  of  the  model  will  be  valid  only  for 

typical  assignments  of  values  for  the  d    ,  a  situation  which  is  familiar 
-^^ n,m' 

in  stochastic  problems.   [Thus  d    =0  (all  n,  m)  is  a  possible  result 
of  a  random  assignment  of  values,  but  it  is  not  typical.   The  probability 
of  this  assignment  vanishes  with  extreme  rapidity  as  N  -^co.^   Instead  of 
restricting  ourselves  to  typical  assignments,  we  could  equlvalently  employ 


a .statistical  distribution  of  assignments  and  make  our  deductions  about 
averages  over  the  distribution.   Such  a  procedure  would  have  some  formal 
advantages,  but  we  feel  that  the  analysis  will  be  clearer  if  we  do  not 
introduce  this  additional  kind  of  average. 

We  wish  now  to  investigate  the  equilibri\;im  thermodynamics  of  our 
model  in  the  classical  limit.   The  virial  expansion  of  the  Helmholtz 
free  energy  per  particle  for  the  true  problem  may  be  written 


P"^  Y^     {a  +  1)'^   p"^   B  .  (2.11) 

ou=l  °' 


Here 


A  =  p"  -|ln|p(2jTfi  P)^'^|  -  ij  is  the  Helmholtz  function  for  free 
particles,  p  is  n/H  ,    B  is  the  Mayer  irreducible  cluster  integral  for 
a  cluster  of  a  +  1  particles,  and  p  =  l/kT,  where  k  is  Boltzmann's  con- 
stant and  T  is  absolute  temperature.   Equation  (2.11)  is  formally  exact 
in  the  limit  il,    N  -^   ™  .   In  the  true  problem  it  makes  no  difference, 
of  course,  which  of  the  N  particles  are  assumed  to  be  in  the  cluster  for 

which  B  is  calculated:  the  interaction  of  all  pairs  is  identical.   If 
a 


that  (2.11)  is  still  valid  provided  that  B  is  reinterpreted  in  the 
following  way:   It  is  the  average  value  of  the  cluster  integral  when 

the  latter  is  calculated  for  all  possible  choices  of  the  a  +  1  particles 

12 
from  among  the  N  particles  in  the  system. 

Let  us  write 

f(x)  =  exp[-pV(x)]  -  1,   f"'"'(x)  =  exp[-pv"''"(x)]  -  1.       (2.12) 


10   - 


Then  the  model   cluster  integral  B,     [Fig.    la~|    is 

B,    =  W"^  n~^     y~]  f  f    f'''"'(x  -   X    )d\  d\   ,  (2.15) 

1  j—J        J  J  ^    n        m^        n        m' 

2 
where  the  integrations  are  over  Q,    and  we  have  replaced  N(N-1)  by  N 

in  anticipation  of  the  limit  N  -^  <»  .   In  the  ladder  model  we  have 

i^'^^ix)    =  f(x  -  d    ). 
Therefore,  since  V(x)  is  cycliC;,  we  find 

\=      I      f(x)d\,  (2.14) 

which  is   identical  with  the   result   for  the  true  problem. 

Let  us   now  assime   that  V(x)    has   a  finite   range   r     and  is   negligible 
for    |x|    >  r    .      (Here  x  is  measured  modulo  displacement  by  a  cyclic  period.) 
The  next  irreducible   cluster  integral    [Fig.   Ibj    is 

n,m,l      '-'  ^  '^ 

X  f^'-(x^-   xJdWx,    .  (2.15) 

Contributions  to  this  integral  can  arise  only  from  points  which  simul- 
taneously satisfy 


Ix  -  x  -  d   I  =  r^;   Ix  -  x,-  d   , I  S  r^,   Ix.-  x  -  d     S  r„. 
'  n   m   n,m'     0    '  m   I        m,  t '     0    '  t   n   i.n' 


0 
(2.16) 
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which  range  over  the  entire  cyclic  volume^  it  vrill  be  impossible  to  satisfy 

(2.16)  for  most  triads  n,  m,  I.      In  fact,  given  a  typical  assignment  of 

the  d's,  it  is  clear  that  (2.l6)  can  be  satisfied  only  for  a  fraction  of 

all  triads  which  is  of  order  r  /fi  .   It  follows  that  Bp  vanishes  as  r     /fl 

in  the  limit  fl   — >  oa,    N  — >  »  ,   Similar  considerations  show  that  any  given 

B   (a  >  1)   also  vanishes  in  the  limit.   The  contribution  to  B  of  each 
a  a 

irreducible  cluster  diagram  with  a  +  1  particles  and  7   links  vanishes  as 
(r  /q)'       .     Actually,  the  condition  we  placed  on  V(x)  is  stronger  than 
needed  to  obtain  this  result.   It  is  sufficient  that  f(x)  be  bounded  and 
that   /  |f(x) |d  X  be  finite  in  the  limit. 

On  the  basis  of  the  preceding  paragraph,  let  us  assume  that  the  total 
contribution  to  (2,11)  from  all  B   (a  >  1)   vanishes  in  the  limit.   This 
is  a  non-trivial  assumption.   It  involves  a  deep-lying  convergence  ques- 
tion wlxLch  we  shall  discuss,  in  its  quantum-mechanical  form,  in  Sec.  5.1, 
and  at  length  in  the  following  paper.   The  essential  point  is  that  the 
number  of  irreducible  diagrams  of  order  a  is  enormous  for  a  ~  0(W).   For 
the  present,  we  shall  simply  adopt  the  assiaiption.   An  equivalent  assump- 
tion will  be  implicit  in  the  discussion  of  the  further  classical  models  of 


3"^a(p)p  ,  (2.17) 


where  -2a(3)  is  the  right  side  of  (2.l4),   The  corresponding  equation  of 
state  is 

p  =  3"^p[l  +  a(p)p]  ,  (2.18) 

where  we  define  the  pressure  by  the  relation  p  =  p''(3A/dp)p. 
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Equations  (2.17)  and  (2.l8)  exhibit  several  properties  of  interest. 
First,  we  note  that  if  V(x)  is  non-negative  everywhere,  then  A  -  A  is 
non-negative  for  all  p  and  goes  to  zero  as  p  —^  «>  „   This  is  consistent 
with  our  previous  finding  that  the  model  potential  energy  is  always  non- 
negative  for  such  V(x).     The  present  result  provides  some  reassursince 
as  to  the  validity  of  our  formal  procedures.   A  second  property  shows 

up  most  clearly  if  we  take  V(x)  to  be  a  hard-sphere  potential  of  range 

2     3 
r  •   Then  we  have  a(p)  =  —  rt  r_^  .   Now  if  we  increase  p  without  limit, 

we  see  that,  in  contrast  to  the  true  problem,  the  ladder  model  exliibits 

Ik 
no  saturation;  the  free  energy  and  pressure  continue  to  rise  smoothly. 

A  tliird  fact  of  interest  is  that  if  V(x)  is  negative  ainywhere,  we 
have  p~  a(p)  — $►  -  <»  as  p  — >  °o  .   Thus  we  have  A  -^  -   «>  for  any  p,  which 
indicates  that  there  is  no  lower  bound  to  the  potential  energy  per 
particle.   Furthermore,  we  have  (5p/Sp)   <  0,  for  any  given  p,  if  the 
temperature  is  low  enough.   This  suggests  that  the  system  then  would  be 
unstable  to  collapse,  and,  since  there  is  no  saturation,  that  the  col- 
lapse would  be  catastrophic  once  it  occurred.     In  the  case  of  poten- 
tials with  an  attractive  part,  the  ladder  model  offers  the  possibility 
of  a  valid  approximation  to  the  true  problem  only  above  a  critical 
temperature  for  each  p.   It  should  be  viewed  with  suspicion  even  above 
this  temperature. 

It  is  clear  from  (2.17)  aJ^d.  (2.l8)  that  the  ladder  model  repre- 
sents an  extremely  rudimentary  approximation  to  the  classical  true 
problem.   The  interest  of  these  results  lies  in  the  fact  that  they 
represent  classical  limits  for  the  fermion  and  boson  ladder  models  which 
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■we  shall  introduce  in  Sec.  3.   These  models  are  non-trivial.  A  similar 
interest  attaches  to  the  further  classical  results  to  be  presented  in 
the  remainder  of  Sec.  2. 

2.3.  Ring  Model 

Instead  of  adopting  (2.9),  let  us  now  specialize 

i  .    =  expri(0  ,  +  G   ,  )I,  (2.19) 

■^n^mjk     -^  ^  n;k    m; -k  r  ^   ^' 

e  ,  =  -0   ,  ,  (2.20) 

n;k     n;-k'  ^    ' 

where  the  9  ,  are  real  phases.   Again  we  see  that  (2.8)  is  satisfied. 
Let  us  give  the  0    values  chosen  at  random  in  the  interval  (0,2jt). 
The  choice  is  to  he  made  independently  for  each  pair  of  indices  n,k, 
subject  only  to  (2.20).  We  shall  call  the  result  the  ring  model,  for 
reasons  which  will  become  clear  shortly.   It  represents  a  rather  more 
drastic  mutilation  of  the  true  Hamiltonian  than  does  the  ladder  model. 
Because  the  phases  0     fluctuate  randomly  as  k  changes,  the  present 
V  '  (x)  are  strange  potentials  which  spread  out  irregularly  over  the 
entire  cyclic  cube. 

The  ring  model  Hamiltonian  also  has  a  boundedness  property  in  common 
with  the  true  Hamiltonian,  but  a  different  one  than  we  noted  for  the 
ladder  model.   Let  us  formally  define  a  self-interaction  potential  by 
extending  (2.6),  (2.7),  and  (2.19),  without  change,  to  the  case  n  =  ra. 
Then  we  may  rewrite  the  model  interaction  Hamiltonian  in  the  form 

H.  =  I  m  V^'"^(x^,-  x^)  -  \   NV(0),  (2.21) 

n,m 
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where  the  summation  now  admits  n  =  m.   Using  (2.6),  (2.19),  and  (2.20), 

we  find 


Z]  exp[i(k.x^^  +  e^.j^)]-  (2.25) 


Now  suppose  that  V  is  non-negative  for  all  k.   Then  ^ |  V  P  P    is  a 

k 
non-negative  operator,  and  it  follows  that  the  expectation  of  the  poten- 
tial energy  per  particle  in  any  quantum-mechanical  state  is  bounded  from 
below  by  -  I  V(0). 

In  the  true  problem  all  the  0  ,  are  zero,  and  P,  is  a  density- 
^  n;k         '  k 

operator  Fourier  component,  as  introduced  by  Pines  and  Bohm   and  others. 
In  the  ring  model,  we  may  call  p   an  effective  density  component. 


exclusive,  but  they  do  not  imply  each  other.   The  bounds  we  have  derived 
therefore  suggest  that  the  ladder  and  ring  models  have  inequivalent 
domains  of  validity.   Consider,  for  example,  the  modified  coulomb  poten- 
tial 

\  =  n'^^rte^  exp(-a|k|)k"^         (lk|  g  l/l), 

(2.2^) 
V^  =  0  (|k|  <  I/O, 

where  a  and  I   provide,  respectively,  a  short-range  and  long-range  cut-off. 
V  is  non-negative  and  V(o)  is  finite.   Thus  the  ring  model  Harailtonian 
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has  a  finite  lower  bound  per  particle  and  may  be  expected  to  yield  healthy 
results.   On  the  other  hand^  V(x)  g  0  is  not  satisfied  for  large  x,  and 
we  cannot  make  a  similar  prediction  for  the  ladder  model. 

If  we  let  I  -^    "> ,    then  V(0)  approaches  a  finite  limit,  and  we  con- 
clude that  the  long-range  character  of  the  coulomb  potential  should  not 
pose  difficulties  for  the  ring  model.   In  the  limit  a  — >  0,  however,  we 
have  -  -p  V(0)  ^  -  <»  .   In  the  true  problem,  the  V  for  very  high  k  give 
a  purely  repulsive  contribution  to  V(x)  and  cannot  actually  cause  H.  to 
be  unbounded  from  below.   In  the  ring  model,  however,  the  very  high  k  give 
rise  to  attractive  as  well  as  repulsive  regions  in  the  T  '    (x),  because 
of  the  fluctuating  phases  of  the  V  '    .   Thus  we  may  anticipate  trouble 
in  the  limit  a  —^  0.   We  shall  see  shortly  that  it  actually  occ\irs,  at 
least  in  the  classical  case. 

The  B  may  be  evaluated  for  the  ring  model  by  expanding  the  f  '  '(x) 
as  power  series  in  -p,  expanding  the  V  '  (x)  in  Fourier  series,  and  then 
performing  the  space  integrations.  We  thereby  find 


,"'■"  + i(-0)"  E  v^"'"v./' 


k,k' 


(2.25) 


By  (2.19)  and  (2.20),  we  find 


n,m  _  n,m^  '^^'^  =  V  V   , 

^0    ~  0'      k    -k      k  -k 


i6 


Thus  the  first  two  terms  on  the  right  side  of  (2.25)  are  iinaffected  by 
the  averaging  over  n  and  m.   All  the  higher  terms ^  however,  involve 
phases  which  fluctuate  randomly  as  n  and  m  are  varied,  except  when  all 
the  summed  indices  k,  k',  ...  are  either  equal  and  opposite  in  pairs  or 
zero.   The  consequence  is  that  none  of  the  higher  terms  makes  a  contri- 
bution to  B  in  the  limit  N  -^  <»,  Q   — >  «>  . 

We  shall  illustrate  by  considering  the  term  containing  (-p)  . 
The  random  phase  of  the  summand  is 

(9  ,  +0  ,)+(9  ,,+9   ,,)+(9   ,  ,,+9  ,  ,,). 
n;k    ra;k    ^  n;k'     m;-k'     ^  n;-k-k'    m;k+k' 

By  (2.20),  this  expression  vanishes  for  k  =  0,  k'  =0,  or  k  +  k '  =0. 
However,  it  follows  readily  from  the  restrictions  on  V(x),  stated  after 
(2.3),  that  the  total  contribution  to  B  from  these  restricted  wave- 
vector  combinations  vanishes  in  the  limit.   For  the  remaining  wave-vector 
combinations,  the  phase  of  the  summand  changes  at  random  with  change  of 
n  and  m.   For  a  given  k  and  k',  the  averaging  over  n  and  m  therefore  re- 
duces the  contribution  to  B,  by  a  factor  '^  W   =  l/  /(N  )   from  its  value 
in  the  true  problem.   The  consequence  is  that  the  total  contribution  of 
the  (-p)   term  vanishes  in  the  limit.   Similar  arguments  show  that  all 
the  higher  terms  vanish  also.   Thus  we  have 


k 


(2.26) 


where  we  note  V  =  V   .   In  obtaining  (2.26),  we  use  the  fact  that  the 
expansion  of  f  '  (x)  in  powers  of  -p  is  absolutely  convergent  for  all  p, 
if  V(x)  obeys  the  restrictions  imposed  after  (2.5). 
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The  higher  B  may  be  evaluated  by  similar  analysis.   The  result 
is  that  the  only  irreducible  Mayer  diagrams  which  give  nonvanishing 
contributions  in  the  limit  are  the  ring  diagrams,  the  first  three  of 
which  are  shown  in  Figs,  lb  and  Ic.   The  surviving  contributions  from 
the  ring  diagrams  give 

B,  =  i  I]  (-P)"-'"  n"  (V^r'  (a.  2).        (2.27) 

The  surviving  contributions  arise  as  the  products  of  the  terms  OC  -p 
in  the  expansions  of  all  the  f  factors  occurring  in  the  ring  diagram 
integrands.   To  see  how  they  survive,  consider  Fig.  lb.   The  surviving 
contribution  from  this  diagram  is 

ir'  e;  e  (-3)'. ■>%,"- v,-'v,'-, 

n,m,  t    k 
and  the  phase  of  the  summand  is 


^  n;k    m;-k    ^  m;k    t;-k    ^  t;k    n; -k  ' 


which  vanishes  by  (2.20). 

Inserting  (2.26)  and  (2.27)  in  (2.11),  and  performing  the  sum  over 


A-Aq  =  \   pHVq  -  I  (3P)"^7"^  X] 


(3Pi;V^_  -  ln(l  +  3PiW^) 


(2.; 


It  is  of  Interest  to  compare  (2.28)  with  the  well-known  result 
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I  pp'^nf^-  |(pp)"^n"^  Y]    [-P"^k  "  ^"^^  "  P"^k^ 


I  P^^Vl'  (2.29) 


)  I  fj^  exp(ik-x); 


17 
which  Montroll  and  Mayer   obtained  by  summing  all  the  ring  diagrams  for 

the  true  problem.   We  shall  find  that  (2.28),  and  not  (2.29),  represents 

the  classical  limit  of  the  quantum-mechanical  ring  summation  to  be  carried 

.  .   ^    1  18 
out  m  Sec.  4. 

From  (2.28),  we  see  that  A-A  is  bounded  from  below  (V  s  o)  by 


Y]   \  =  -^v(o). 


which  agrees  with  the  rigorous  bound  we  have  previously  found  for  the 
ring-model  potential  energy.   As  in  the  case  of  the  ladder  model,  this 
provides  some  reassurance  as  to  the  validity  of  our  formal  procedure. 
It  is  clear  that  if  V  =  0,  then  A-A  will  actually  approach  the  absolute 
lower  bound  as  3  — >  oo.   [The  In  term  in  (2.28)  gives  a  vanishing  contri- 
bution in  this  limit.] 

If  V  is  given  by  (2.2^),  we  find  that  A-A  converges  in  the  limit 
I  -~>  <xi   and/or  p  — >«.   This  supports  our  anticipation  that  the  long-range 
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part  of  the  potential  should  not  pose  difficulties  for  the  ring  model. 
It  sho\ild  be  noted  that  the  derivation  of  (2.28)  with  potential  (2.2'+) 
requires  that  t   be  kept  finite  until  after  the  limit  Q,   — >  «  is  taken. 
Otherwise,  the  assumption  V,  =  0(n   )  is  violated  for  very  low  k. 

If  now  we  take  a  =  0,  we  find  A-A  — >  -  «  as  3  ^  «>.   Thus,  as  ve 
anticipated  might  be  so,  the  ring  model  is  not  an  admissible  approxima- 
tion in  this  case.   The  situation  m.ay  be  substantially  improved  in  the 
quantum  theory,  however. 


2.k     Random- Coupling  and  Hartree-Fock  Models 

We  wish  now  to  examine  two  simpler  distinguishable-particle  models, 

In  common  with  the  ladder  and  ring  models,  they  are  of  interest  because 

their  classical  thermodynamics  represents  limits  for  corresponding 

fermion  and  boson  models. 

Let  us  now  specialize  ^         ,  to  the  form 
"^n^mjk 


i         .    =   exp(i0    ,  ),  (2.50) 

n,m;k     ^^      njffljk'^'  ^   -^  ' 


with 


G    ,  =  -0    ,  ,     9    ,  =  -9     ,  ,  (2.51) 

n,m;k     m,n;k        n,m;k     n,m;-k  ^     ' 


where  the  9      are  real  phases.   Again  (2.8)  is  satisfied.   Let  us 

give  the  9    ,  values  in  the  interval  (O,  2n)   determined  by  indenendent, 
^        n,m:k  \    f        /  j  ^  > 

random  choices  for  all  the  combinations  of  indices, subject  only  to  (2.51). 
We  shall  call  this  the  random- coupling  model. 
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There  appears  to  be  no  lower  boimd  to  the  potential  energy  per 
particle  in  the  random  coupling  model  if  we  take  the  limit  W  ^  oo. 
Consequently;  we  cannot  be  sure  that  the  model  has  any  thermodynamic 
validity.   We  shall  return  to  this  question  in  a  moment. 

The  B  for  the  random- coupling  model  may  be  determined  by  the  same 
formal  procedure  as  we  used  for  the  ring  model.   The  result  is  that  B 
has  the  value  (2.26)  and  that  all  the  higher  B  vanish  in  the  limit 


The  results  for  A  and  p  are 


A-Aq  =  i  PQVq  -  J  ppn  XH  \^  (2-52) 

p  =  p'^P  +  I  p^%  -  J  PP^^  Z]  \^-  (2-33) 

k 


If  (2.17)  is  expanded  in  powers  of  -fi,    and  then  expressed  in  terms  of 
Fourier  coefficients^  it  is  easy  to  verify  that  (2.52)  represents  that 
part  of  (2.17)  which  is  also  contained  in  (2.28).   That  is  to  say,  the 
only  contributions  to  A  which  survive  in  the  random- coupling  model  are 
those  which  survive  in  both  the  ladder  and  ring  models. 

From  (2.53)  we  see  that  if  p  and  (3  are  high  enough  we  have 
(Sp/Sp)   <  0,  regardless  of  the  form  of  V(x).   The  instability  to  col- 
lapse thereby  indicated,    is  associated  with  the  lack  of  a  lower  bound 
to  the  model  potential  energy.   However,  if  p  and  p  are  low  enough, 
(2.55)  suggests  that  the  random- coupling  model  may  have  a  stable  thermo- 
dynamics.  We  choose  to  regard  that  indication  with  caution. 

Our  final  model  is  the  Hartree-Fock  model,  which  we  construct  by 
taking 
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V/'^^^O         (k^O).  (2.54) 


This  is  a  zeroth  model  in  the  sense  that  there  are  no  randomly  chosen 
parameters  at  all.   It  corresponds  simply  to  having  each  particle  move 
in  the  \miform  potential  obtained  by  averaging  the  true  fields  of  the  other 
particles  over  all  possible  configurations.   Only  B^  is  nonvanishing  in 
the  limit  Q.   — >  «>,  and  we  have 

A-Aq  =  I  pAVq.  (2.55) 


5.   MODELS  FOR  FERMIONS  AND  BOSONS 
5.1.  Nature  of  the  Models 

The  models  described  in  Sec.  2  involve  interaction  potentials  which 
are  different  for  different  pairs  of  particles.   They  are  therefore  mean- 
ingless for  indistinguishable  particles.   In  order  to  construct  stochastic 
models  for  fermion  and  boson  systems,  let  us  replace  (2.1)  by  the  second- 
quantized  true  Hamiltonian 


^0=  5  \\\'  ^5.1) 

k 

«i  =  i  fli   \-s\+P,r-^s\Ws•  (5-2) 

kprs 


Here  q   and  q  are  fermion  or  boson  creation  and  destruction  operators 

1  2 

for  momentum  k,  e  is  the  free-particle  energy  -^.'k.   ,    and  we  take  n  =  1. 
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The  commutation  relations  are 


vhere  the  plus  sign  is  for  fermions  and  the  minus  for  bosons. 
As  the  general  model  interaction  Hamiltonian;,  we  take 

«i  =  i     XJ    \-s^kprs\+p,r+s\\'Vs    '  ^^A) 

kprS  Jr-        Jr-; 


where  the  jz^     are  c-number  parameters  which  play  a  role  analogous 
to  that  of  the  (i  ,  .   We  leave  H^  unaltered.   In  correspondence  to 

(2.8),  we  impose  the  conditions 

^kprs  =  ^srpk*  '     «^kprs  =  ^pksr  '  ^5-5) 

The  first  of  these  relations  insures  the  Hermiticity  of  H. .   The 
second  is  suggested  by  the  invariance  of  (5-2)  to  the  'particle 
exchange'  (k,s)  - — ^  (p,r). 

We  shall  obtain  the  fermion  and  boson  versions  of  the  ladder,  ring, 
random- coupling,  and  Hartree-Fock  models  by  making  specialized  stochas- 
tic assignments  of  values  to  the  jz^^  ^^ .   The  <i>^  will  all  have  unit 
modulus  in  the  models  we  shall  examine  here,  except  in  the  Hartree-Fock 
model,  where  most  of  the  (j)        ^   will  vanish.   In  Sees,  h   and  5  we  shall 
develop  an  appropriate  propagator  formalism  for  the  fermion  and  boson 
models  and  find  closed  equations  which  determine  the  propagators  for 
each  model.   We  shall  put  off  until  Sec.  6  a  demonstration  of  the  rela- 
tions among  the  fermion,  boson,  and  distinguishable-particle  models. 
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5.2.  Ladder  Model 

To  construct  the  fermion  or  boson  ladder  model,  we  take 

with 

e,  =  e  ,  .  (5-7) 

kp    pk  \ ^  I / 

We  then  determine  the  real  phase  0   for  each  pair  k,p  by  a  random 
choice  in  the  interval  (O,  2rt).   The  choices  are  all  independent, 
subject  only  to  (5.7). 

Let  us  define  the  quantities 

X(x',x)  =  n~^     2_i  %%   expri(r-x'  +  s-x  +  ©^3)], 


(5.8) 


X^ix\x) 


Q   ^  2_J  %^%^   expr-i(r-x'  +  s-x  +  0^^)1 


In  the  true  problem  (all  9  ^  =  o),  x(x',x)  is  simply  the  two-particle 
amplitude  il/(x' )\|/(x) ,  where  i|/(x)  is  the  destruction  field  in  x  space. 
We  may  call  X(x',x)  the  effective  two-particle  amplitude  in  the  ladder 
model.   By  straightforward  Fourier  analysis,  we  find 

H.  =  I   /"r  V(x-x')x"^(x',x)x(x',x)d\d\'.  (5.9) 

If  V(x)  g  .0  for  all  x,  H.  is  a  positive-definite  operator,  and  it  follows 
that  the  eigenvalues  of  H  are  all  non-negative.   This  is  the  same  boiind 
which  we  obtained  in  Sec.  2  for  the  distinguishable-particle  ladder  model. 


2h 


It  is  not  obvious  that  the  present  model  will  be  an  admissible 
approximation  to  the  true  problem  if  V(x)  is  a  hard-sphere  potential, 
although  this  was  clearly  the  case  for  the  distinguishable-particle 
ladder  model.   A  sufficient  condition  for  admissibility  would  appear 
to  be  that  the  relation 

X(xSx)|$>  =  0,  (|x-x'|  %    T^),  (5.10) 

fi 

fold  of  states  ^  in  the  model  as  in  the  true  problem.   We  have  not  in- 
vestigated this  question. 

3-5.  Ring  Model 

To   construct  the   fermion  or  boson  ring  model,   we  take 

^,  =   expli(a       +9      )1,  (3.11) 

"^kprs  -^1    ^   ks  pr  J 


(3.12) 

We  fix  the  real  phases  0   by  random  choices  in  the  interval  (O,  2jt ) , 
subject  only  to  (3.12) . 

In  analogy  to  (2.23),  let  us  introduce  the  effective  density- 
component  operators 

P  =  /    q/q,    exp(i0,  .   ).  (3.13) 

m    4r—  ^k  ^k+m   ^^   k,k+m' 
k 

By  (3-12),  they  satisfy  p  =  p   ^.   It  follows  from  (3-3)  and  (5-12)  that 
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for  either  fermlons  or  bosons  the  ring  model  H.  may  be  rewritten 


where 


As  in  (2.22),  we  note  that  the  first  terra  on  the  right  side  of  (3.1^) 
is  a  positive-definite  operator  if  V  S  0  for  all  k.   It  follows  that 
the  fermion  and  boson  ring  models  exhibit  the  same  lower  bound  on  H 
as  did  the  distinguishable-particle  ring  model. 

A  third  model,  which  has  no  analog  in  the  distinguishable  case, 
may  be  constructed  by  taking 


Ws  =  'M^^\r  ^  %s^  (5.15) 


and  requiring  the  9   to  obey  (3.12).   We  may  call  this  the  exchange 
model.   Hybrid  models  may  also  be  constructed,  by  taking  ^,     as  a 
linear  combination  of  the  forms  for  the  ladder,  ring,  and  exchange 
models.  We  shall  not  discuss  these  cases  in  this  paper. 


3.^.  Random- Coupling  and  Hartree-Fock  Models 

To  construct  the  fermion  and  boson  versions  of  the  random- couplir 
model,  we  take 


W  =   ^^P^^Vs^'  ^5.16) 
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e,    =  -9   ,  ,  0,    =  0  ,   ,9,    =  9  ,   ,         (3-17) 
kprs     srpk    kprs    pksr    kprs    pkrs 


and  fix  the  phases  9     by  independent  random  choices  for  each  com- 
bination of  indices  (k,p,r,s),  subject  only  to  (J.IT).   As  in  the 
distinguishable-particle  random  coupling  model;  there  appears  to  be  no 
lower  bound  on  the  potential  energy  per  particle  when  the  system  is 
infinite.   This  suggests  that  results  obtained  from  the  random- coupling 
model  be  viewed  with  skepticism.   It  perhaps  should  be  pointed  out  that 
our  random- coupling  models  are  unrelated  to  the  random-phase  approxima- 
tion employed  by  Pines  and  Bohm   and  others.   We  make  no  assimiption 
about  the  phase  correlations  among  the  dynamic  variables. 

The  fermion  or  boson  Hartree-Fock  model  is  given  by  the  assignment 

<!>       ,     =  i,    ,       =   \,      I  =0     (k  /  r  or  s).  (5.18) 

'^kppk   '^kpkp    '   ^kprs 

This  model  contains  no  randomly  chosen  parameters.   We  shall  see  in 
Sec,  5  that  it  yields  simply  the  Hartree-Fock  approximation  to  the  true 
problem,  in  the  limit  Q  — >  ». 

h.         TEMPERATURE-DOMAIN  PROPAGATOR  FORMALISM 

The  equilibrium  statistical  mechanics  of  the  fermion  and  boson 

models  can  be  investigated  most  neatly  by  means  of  the  temperature- 

2,5,6,8 
domain  propagator  formalism  and  its  associated  diagram  technique. 

Since  our  models  differ  from  the  true  Hamiltonian  only  by  the  replacement 

V   6        ->  V   (z5,    6,       ,  C+'l) 

k-s  k+p,r+s     k-s'^kprs  k+p,  r+s' 
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the  existing  techniques  may  be  taken  over  with  only  minor  and  obvious 

changes.  We  shall  summarize  the  resulting  formalism  in  the  present 

19 
section. 

Let  us   define  the   temperature-domain  propagator  S,  (u,u')   by 
Sj^(u,u')    =    -     {T[qj^(u)q^'^(u')]^  (u,u'    real),  {k.2) 

where 

/    X  uH        -uH  if    .  uH     t   -uH  ,,     ,^ 

\{^)   =   ^     \^  '  \   (u)    =   e      q^^  e         .  (4.5) 

The  ordering  operator  T  is  defined  by 

Trqj^(u)qj^^(u')j  =  qj^(u)q^'^(u' )  (u  >u'). 


(4.M 


Tqi^(u)qj^  (u')l  =  +  \   (u')qj^(u)        (u  ?  u') 


In  (4.'+)  and  in  all  subsequent  expressions  where  a  plus-minus  or  minus- 
plus  sign  occurs,  the  upper  sign  refers  to  fermions  and  the  lower  to 
bosons.   The  brackets  \^  denote  an  average  over  the  grand  canonical 
ensemble.   For  any  operator  B, 

<B>  .  Tr{e-P^M)3|   /  ^^  {e-3(H-MN)  |  ^  ^^^^^ 

where  N  is  defined  below  (J.l^)  and  |a  is  the  chemical  potential. 
The  propagator  has  a  Fourier  expansion  of  the  form 

Sj^(u,u')  =  p-^    5Z   Sj^(^^)expR  (u'-u+6)1,  (i+.6) 


^^    =   ^   +  ijt(2a  +  l)p  (fermions), 


^     =  H  +  2ijta|3  (bosons)^ 


c+.t: 


and  a  takes  all  integer  values.   The  quantity  5  is  an  infinitesimal  real, 
positive  number.   We  shall  call  S  (^  )  a  propagator  also,  and  we  shall 
call  C,^    an   'energy'. 

The  complete  thermodynamic  behavior  of  the  system  can  be  obtained 
from  S  (^^).   The  mean  number  of  particles  at  a  given  temperature  and 
chemical  potential  is 


k  a 

where  the  W   =  <f^,/       are  the  mean  occupation  numbers.  The  mean  energy 

also  has  a  direct  expression.   Using  (3-1);.  (50)j>  {3-^) ,  and  (^.J)?  we 
find 

±  E  K(u,u.)/Su.]^,^^=  <H„>  .2    (k.)     ,  (K.9) 

and  we  note  that  <^hA  =  -     Z_j  ^y.\i^,^)  ■      Then  it  follows  from  {h.6) 

20 
that  the  mean  energy  E(0,h)  is  given  by 


E(p,^i)  =  ±1  P'^  YZ     (V  ^a^^^^a^^^P^^a^)-  ^^"^^^ 

k,  a 

The  entropy,  pressure,  and  other  thermodynamic  quantities  can  be 
found  from  N(p,ij)  and  E(p,n).   [Alternatively,  the  thermodynamic  poten- 
tial may  be  obtained  from  S,  ( ^  )  by  an  integration  over  an  interaction 
strength  parameter.  J 
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The  propagator  for  free  particles  is 


\'°'«,)  =  U,  -  e^)-\  (k.Xl) 


The  coupled-particle  propagator  S,  (^  )  may  be  expressed  in  terms  of 


(0), 

following  rules: 


S^   (^  )  hy  a  linked-diagram  expansion  constructed  according  to  the 


1.  Call  the  diagram  part  shown  in  Fig.  2a  a  vertex.   A  vertex  con- 
sists of  two  solid-line  junctions  connected  by  a  dashed  line.   'Line' 
hereafter  will  mean  solid  line  except  where  noted. 

2.  For  each  positive  integer  n,  take  n  vertices  and  join  incoming 
with  outgoing  lines  in  pairs  so  as  to  form,  just  once,  all  possible 
distinct,  linked  diagrams  with  just  one  incoming  and  one  outgoing  ex- 
ternal line.   Linked  diagrams  are  those  which  do  not  consist  of  dis- 
connected parts.   External  lines  are  those  which  leave  or  enter  the 
diagram.   In  reckoning  distinctness,  the  n  original  vertices  are 
considered  indistinguishable  and  the  two  solid-line  junctions  in  any 
vertex  are  considered  indistinguishable.   However,  incoming  lines  are 

distinct  from  outgoing  lines.   Examples:  Figs.  6a  and  8a  are  distinct, 

21 
but  Figs.  6a  smd  9a  are  not. 

5.   Label  the  external  lines  i,d.th  momentum  k.   Label  the  internal  lines 
with  momenta  k '  ,  k",  ...  .   In  addition,  associate  the  'energy'  t,     with 
the  external  lines,  and  'energies'  ^  ,,  ^^„  ,  ...   with  the  internal 
lines . 
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k.        With  each  llne^  external  or  internal^  associate  a  factor  S    ^^y.'^' 
where  p  is  the  momentum  labeling  the  line  and  t  is  the  associated 
'energy'.   Special  case:   Include  an  additional  factor  exp(C5)  if  the 
beginning  and  termination  of  the  line  are  in  the  same  vertex. 

5.  With  each  vertex,  at  which  are  joined  lines  with  momentum  labels 

p,  q,  T,    s  and  respective  'energies'  L,    ^  ,  t,^,    ^^   as  shown  in  Fig.  2b, 
associate  a  factor 

-B'"'"  2$     V    5        5^    J 
^   '^pqrs   p-s   p+q,  r+s  b+c,d+e 

Note:   This  rule  is  ambiguous  with  respect  to  the  exchange  (p, s)  ^r±  (q, r) . 
By  (2.5)  and  (5-5),  however,  the  factor  called  for  by  the  rule  is  in- 
variant in  value  under  the  exchange. 

I 

6.  For  fermions  only:   Associate  with  each  diagram  a  factor  (-1)  ,  where 

I   is  the  number  of  closed  solid-line  loops  in  the  diagram. 

7.  To  form  the  contribution  of  a  given  diagram,  multiply  together  all 
the  factors  introduced  by  rules  k,    5,  and  6  and  then  sxaa   over  all  the 
momenta  k',  k",  ...   and  'energies'  ^  ,  ^  ^  n  ^  • • •   associated  with  the 

internal  lines . 


8.   To  form  S,  (^  ),  first  sum  the  contributions  given  by  rule  7  for  all 
the  diagrams  admitted  by  rule  2.   Then  add  the  contribution  S^^   ( ^a^ ' 
which  is  associated  with  the  zeroth  order  (n  =  O)  diagram  Fig.  Ja. 


We  shall  call  the  expansion  for  S  (^  )  generated  by  rules  1-8 
the  primitive  linked-  diagram  expansion.   A  more  compact  expansion,  which 
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we  shall  call  the  irreducible  linked- diagram  expansion,  may  be  constructed 
by  replacing  rules  2  and  h   with  the  following  altered  rules: 

2' .   Retain  only  those  diagrams  admitted  by  rule  2  which  do  not  contain 
self-energy  parts.  We  shall  call  these  irreducible  diagrams.   (  A  self- 
energy  part  is  a  part  of  a  diagram  which  contains  at  least  one  but  not 
all  the  vertices  and  which  is  connected  to  the  rest  of  the  diagram  by 
jiist  one  outgoing  and  one  incoming  line . ) 

k' .      In  rule  k,    replace  each  factor  S  ^   (L)  ^y  a  factor  S  H^)- 
Exception:   With  the  outgoing  external  line  associate  the  factor  S    (^  ) 
as  before. 

The  primitive  diagram  expansion  gives  S,  (^  )  as  an  infinite  s\am  of 
integrals  (we  consider  the  case  Q,  -^  oo)  over  the  known  quantities 
S^   (^  )•   In  contrast,  the  irreducible  expansion  is  really  an  infinite- 
series  integral  equation  for  S  (^  ).  A  convenient  way  to  express  the 
irreducible  expansion  is  as  follows:   We  define  H  (^  )  by 

Then  the  irreducible  linked-diagram  expansion  for  M.  (^  )  is  the  same 
as  that  for  S,  (^  )  ^except  for  the  changes  expressed  by  the  following 
further  rule  alterations: 

h" .      In  rule  k' ,    omit  entirely  the  factors  for  the  two  external  lines. 

8' .   In  rule  8,  omit  the  contribution  of  the  zeroth-order  diagram. 
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The  use  of  the  rules  will  "be  illustrated  by  the  examples  treated 

in  Sec.  5.   If  ve  talce  iz^,     =1  for  all  k^  p,  t,    s,    then  the  rules  we 
-^  "^kprs 

have  given  yield  the  established  propagator  expansions  for  the  true 
problem. 


5.    PROPAGATOR  EQUATIONS  FOR  THE  MODELS 

5.1.  Underlying  Assumptions 

We  shall  now  show  that  our  fermlon  and  boson  models  lead  to  closed 
integral  equations  for  S  ( ^^  )  in  the  limit  H  ^^  «  with  fixed  |-i  or  P.  We 
wish  first  to  state  clearly  the  assumptions  which  underlie  the  analysis: 

1.  As  in  Sec.  2,  we  restrict  V(x)  to  be  a  smooth,  bounded  function 
such  that  |V  I  S  0(k"  ),  k  ->  co. 

2.  We  assume  that  all  the  sums  over  intermediate  momenta  which 
occur  in  any  given  order  of  the  primitive  linked-diagram  expansion  con- 
verge at  infinity.   More  precisely,  we  assume  that,  for  given  k  and  a, 
the  contributions  to  S  (^  )  which  involve  intermediate  momenta  higher 

than  some  given  momentum  k    vanish  as  k    — >  co.   Moreover,  we  assume 
^  max  max 

that  they  vanish  in  a  manner  independent  of  H  as  H  ->  <».   We  believe 
that  this  assumption  actually  follows  from  assumption  1,  but  we  shall 
not  tr>'  to  prove  this  here. 

5.   As  in  Sec.  2,  we  assume  that  V  =  0(n~  )  for  all  k  as  fi  ->  00. 
This  means  that  for  long-range  potentials  a  cut-off  length  I   must  be 
employed  as  in  {2.2h) .      We  take  the  limit  I  ->  «>  only  after  the  limit 
n   — >  00  has  been  performed. 


55 


(0), 

independent  of  il  as  n  ->  oo  with  P  constant .   This  excludes  from  our  pre- 
sent considerations  boson  systems  below  the  Einstein-Bose  condensation 
temperature , 

5.   The  final  ass\;uiiption  involves  a  deep-lying  convergence  question 
corresponding  to  that  which  arose  in  Sec.  2.  We  shall  see,  for  each  of 
the  models,  that  large  classes  of  diagrams  give  a  vanishing  contribution 
to  S,  (^  )  in  the  limit  Q,  —>  <^,    up  to  any  given  order  of  diagram.   On 
this  basis,  we  shall  assume  that  these  classes  do  not  contribute  when 
simmed  to  all  orders.   It  will  not  follow  from  oior  analysis  that  this  is 
actually  so.   The  reason  is  that  the  diagrams  of  order  N  {=(£1)    and  high- 
er are  enormous  in  number  for  large  Q,  and  it  will  not  be  clear  that 

cancellations  due  to  the  randomness  of  the  i,  will  suppress  the 

"^kprs        -^-^ 

contribution  of  these  diagrams  as  they  do  contributions  of  finite  order. 

The  necessity  for  assumptions  5  and  h-   will  be  eliminated  by  the 
generalized  treatment  given  in  the  following  paper.   There  we  construct 
models  for  systems  of  any  size  and  obtain  closed  propagator  equations 
without  taking  the  limit  Q  — 5»-oo,   The  final  equations  are  identical  with 
those  to  be  derived  here,  and  they  justify  the  latter  in  cases  where 
assumptions  5  and  h   are  not  satisfied.   In  particular,  they  establish 
the  closed  model  equations  for  condensed  boson  systems.   The  generalized 
treatment  provides  much  neater  derivations  of  all  the  results  to  be 
obtained  in  Sec.  5.  We  do  not  employ  it  at  the  outset  because  it  re- 
quires an  unfamiliar  and  more  elaborate  formalism. 
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Assumption  5  also  is  best  examined  by  the  methods  of  the  following 
paper.   The  approach  to  tliis  question  adopted  there  is  to  consider  the 
propagators  in  the  real-time  domain,  rather  than  the  temperature  domain, 
and  to  regard  them  as  the  limits  of  more  general  quantities  (correlation 
and  Green's  functions)  which  are  defined  for  non-equilibrium  as  well  as 
equilibrium.   Linked  diagram  expansions  exist  which  give  the  evolution 
of  the  correlation  and  Green's  functions  forward  in  time  from  a  given 
initial  statistical,  state.   As  in  the  present  case,  these  expansions 
can  be  formally  summed  and  expressed  by  closed  equations  for  each  of 
our  models . 

A  new  feature,  however,  is  that  the  evolution  in  time  can  also  be 
examined  by  an  alternate  method  which  seems  genuinely  independent  of 
perturbation  expansions.   One  can  replace  the  exact  differential  equa- 
tions of  evolution  by  corresponding  difference  equations  involving  small 
time  increments.   In  contrast  to  a  perturbation  expansion  (which  is  akin 
to  a  Taylor  series)  such  a  procedure  should  converge,  as  the  increment 
size  is  decreased  to  zero,  whenever  the  differential  equations  themselves 
are  meaningful.   This  permits  an  examination  of  assumption  5  from  a  new 
point  of  view.   Although  the  analysis  we  shall  present  in  the  following 
paper  is  not  rigorous,  we  feel  that  it  provides  substantial  support  for 
the  validity  of  our  closed  model  equations. 

5.2.   Hartree-Fock  and  Random- Coupling  Models 

The  derivation  of  closed  propagator  equations  is  simplest  for  the 
Hartree-Fock  model.   Consider  the  primitive  linked-diagram  expansion  for 
S,  (^  ).   There  are  two  distinct  first-order  diagrams,  and  they  are  shown 
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in  Figs.  5b  and  3g.  By  the  rules  of  Sec.  h,    the  contributions  of  these 
diagrams  involve  the  factors  ^  and  ^    ,    respectively^  but  no  other 

^  factors.   Therefore,  by  (5.I8),  their  contrihutions  are  identical  in 
the  Hartree-Fock  model  and  in  the  true  problem. 

Hovever,  consider  Fig.  ^a  .   The  contribution  of  this  second-order 
diagram  to  the  primitive  expansion  for  S,  (^  )  is 


*   P"'  g  K^rs*sr^A-s\-V.    '^    %!°\i^)%^°'' U^^ 


XS^(°>(C,)s;°>(c,)sJ°)(5,),      (5.1) 


where  r  =  k  +  p-s,  d=:a  +  b-c,  by  momentum  and  'energy'  conserva- 
tion. By  (5. 18),  the  only  surviving  terms  in  the  sum  are  those  for 
which  either  s  =  k,  r  =  p  or  r  =  k,  s  =  p.  Thus,  there  is  only  one 
free  intermediate  momentiMi,  which  we  may  take  as  p.  Kow  as  Q,  — >  00, 
the  number  of  allowed  momenta  p  in  any  given  volLmie  of  momentum  space 
is  ccr  f2.  It  then  follows  from  assumptions  2,  5,  and  ^1-  of  Sec.  5.I 
that  the  contribution  (5-1)  vanishes  as  D, 

If  similar  considerations  ares^jplied  to  the  rest  of  the  primitive 
linked  diagrams,  it  is  possible  to  verify  the  following  result:   The 
only  diagrams  which  survive  in  the  limit  fi  — >  <»  are  those  for  which 
momentijm  conservation  alone  assures  that  the  (^   product  given  by  the 
rules  of  Sec.  k   consists  wholly  of  factors  of  the  form  i^^     of  ''^pqpq  * 
Such  diagrams  give  the  same  contribution  as  they  do  in  the  true  problem. 
For  every  other  diagram  of  finite  order,  (5.18)  so  restricts  the  summa- 
tions over  intermediate  momenta  that  the  result  vanishes  as  some  positive 


-  36 


It  is  easy  to  see  from  the  rules  of  Sec.  k   that  the  higher  surviving 
primitive  diagrams^  which  we  have  just  specified,  are  simply  those  which 
can  be  obtained  from  Figs.  Jb  and  yc   by  repeatedly  inserting  these  first- 
order  diagrams  into  themselves  and  into  each  other  as  self-energy  parts. 
(See  Fig.  5.)   Now  let  us  invoke  assumption  5  of  Sec.  5.1.      It  then 
follows  that  Figs.  Jb  and  Jc  are  the  only  diagrams  which  survive  in  the 
irreducible  expansion  for  S  (^  )  in  the  limit  H  ^  «>.   We  may  see  this 
by  noting  that  the  primitive  expansion  can  be  recovered  from  the  irre- 
ducible expansion  by  talcing  every  factor  S  which  occurs  in  the  latter  and 
replacing  it  by  its  own  primitive  expansion.   If  diagrams  other  than 
Figs.  5b  and  5c  occurred  in  the  irreducible  expansion,  then  it  is  clear 
that  the  primitive  expansion  for  S  ( ^^  )  thus  recovered  would  contain 
diagrams  other  than  those  we  have  specified. 

Using  the  rules  of  Sec.  ^  to  evaluate  the  contribution  to  H  (^  ) 
of  the  two  surviving  irreducible  diagrams,  we  find 

V^a)  =  P"'  ^     ^±  V  \-p)y^b)^-^^^bS^'  (5.2) 

pb 

where  we  have  noted  (2.3).   The  first  tenn  on  the  right  side  of  (5'2) 
arises  from  Fig.  3b  and  the  second  from  Fig.  5c.   Upon  inserting  (5-2) 
into  (i)-.12),  we  obtain  a  closed  integral  equation  for  S  (^  )  for  the 
Hartree-Fock  model.   Equation  (5.2)  may  be  rewritten  in  the  form 

VU=  C  (V^^  v^,_p)^p^  (5.3) 

where  we  use  (^.8).   The  parts  of  this  result  involving  V  and  Vj^_   are, 


37 


respectively,  the  direct  and  exchange  parts  of  the  effective  potential 
which  is  obtained  in  the  usual  Hartree-Fock  approximation  to  the  true 
problem. 

Let  us  consider  next  the  random- coupling  model.   By  (3.l6)  and 
(5.17).  ve  find 


k      ^    =  ^>      k    ^       =  ^>      k         ^       ,    =  1,    k         jZ$  ,   =  1.      (5. if) 
kppk       '^kpkp       '^kprs^si^k    '    '^kprs'^srkp  ^^   ' 


It  then  follows  from  the  rules  of  Sec.  k-   that  Figs.  5b,  5c,  ^a  and  'i-b 
all  survive  in  the  primitive  expansion  for  S,  (^^)  and  give  the  same  con- 
tributions as  in  the  true  problem. 

However,  consider  Fig.  6a.   The  contribution  which  it  makes  to  the 
primitive  expansion  for  S  (^  )  is  of  the  form 


+  (3     _/ I   kprs'^srr's'  s'r'pk  k+p,  r+s  r+s,r'+s'  r'+s',p+k 

prsr's' 


}     I  (product  of  S^°^  factors).    (5.5) 

(energies) 


There  is  an  identity  among  the  conservation  conditions  given  by  the 
three  Kronecker  symbols.   Consequently  there  are  three  independent 
intermediate  momenta,  which  we  may  take  as  p,  s,  and  s'.   For  special 
values  of  these  momenta,  the  !^   product  is  \anity  by  (5.^).   However,  in 
correspondence  to  the  result  noted  above  for  the  Hartree-Fock  model, 
it  is  easy  to  see  that  these  special  values  give  a  vanishing  contribu- 
tion to  S  (^  )  when  fi  ^  ».   Except  for  these  restricted  momentxim  values. 


-  38  - 

the  <j:>   product  will  have  a  phase  which,  by  (3.l6)  et  seq.,  fluctuates 
at  random  with  change  of  P;,  s,  and  s'  . 

Now  let  us  divide  the  momentum  space  into  small  regions,  of  'volume' 
A  ,  such  that  S    x'^iS)    exhibits  negligible  change  if  p  varies  within  a 
given  region.  By  assumption  4  of  Sec.  5-lj  this  should  be  possible  in 
the  limit  Q,   — >  ».   Now  as  f2  -^^,    the  density  of  allowed  modes  in  momen- 
tum space  is  ^    Q.      Hence,  in  the  summation  over  momenta  in  (5.5),  the 
number  of  terms  for  p,  s,  and  s'  within  given  small  regions  will  be 
^^  (nA)  .   Let  us  consider  the  contribution  from  those  terms  in  which  the 
phase  of  the  ^   product  fluctuates  at  random  with  change  of  p,  s,  and  s'. 
As  we  have  just  noted,  these  constitute  all  but  a  negligible  fraction  of 

the  terms.   Because  of  the  phase-fluctuation,  the  total  contribution 

-3/2 
of  these  terms  will  be  down  by  a  factor  ^ (QA)      '      from  the  value  it 

would  have  in  the  time  problem.   Since  this  factor  vanishes  as  Q  -^  00^ 

it  follows  from  assumptions  2  and  3  of  Sec.  5.I  that  the  contribution  to 

S  (^  )  from  Fig.  6a  vanishes  in  the  limit.   Our  argument,  of  course, 

assumes  a  typical  assignment  of  values  to  the  randomly  chosen  parameters 

ft      (cf.  the  discussion  in  Sec.  2.2). 
kprs 

Similar  analysis  may  be  applied  to  the  Mgher  diagrams  in  the 
primitive  expansion  for  S  (^  ).   The  result  is  that  the  only  diagrams 
which  survive  are  those  for  which  momentum  conservation  alone  assures 
that  the  (j)   product  consists  entirely  of  factors  and/or  factor-pairs  of 
the  forms  shown  in  {'j.h).      Any  other  diagram  gives  a  contribution  which 
vanishes  as  some  negative  power  of  fi  as  Q  -^  ».   The  surviving  diagrams 
can  be  seen  to  be  those  which  can  be  constructed  from  Figs.  Jb,  3c,  ^a. 
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and  h-h  by  inserting  these  same   foxir  diagrams   repeatedly  as   self-energy 
parts,    in  correspondence  to  the   sitiiation  for  the  Hartree-Fock  model. 
It   follows   from  this  that  the   only  siirvlving  diagrams   in  the  irreducible 
expansion  for  S,  (^    )   are  Figs.    5b,    5c,    ha.,    and  ^b.      Again,    assimption  5 
of  Sec.    5.1  is  implicit  in  the   argument. 

Combining  the   contributions   of  the   four  surviving  irreducible  dia- 
grams  according  to  the   rules   of  Sec.   h,   we   find 


\aj  =    I]    (Vo  +    V^_p)Np 


^  P"'     II!    V^.3(V^_3  -^  V        )S   (^^)S^(5^)S  (^^^^_^),  (5.6) 

psbc 


where,  again,  we  note  (^.8). 

If  the  third  of  constraints  (5.17)  were  relajced  in  assigning  random- 
values  to  the  a    ,  we  would  obtain  a  version  of  the  random- coupling 
kprs' 

model  in  which  the  exchange  diagram  Fig.  ^b  did  not  survive.   The  corre- 
sponding expression  for  ML  (^  )  would  not  include  the  term  in  (5'6)  which 

involves  V 

p-s 

5.5.   Ring  and  Ladder  Models 

The  results  described  above  for  the  Hartree-Fock  and  random- coupling 
models  may  be  sxjmmarized  very  simply:   The  only  diagrams  which  survive  in 
the  irreducible  expansion  for  R  (^  )  are  those  for  which  the  associated  ^ 
product  has  the  value  one  for  all  values  of  the  intermediate  momenta 
allowed  by  momentiom  conservation.   The  contributions  of  the  surviving 


ko 


diagrams  have  precisely  the  same  form  in  the  models  as  in  the  true  prob- 
lem.  This  general  result  is  also  true  for  the  ring  and  ladder  models 
imder  the  assumptions  of  Sec.  5'1>   Analysis  similar  to  that  we  have 
described  shows  that  the  irreducible  diagrams  in  which  the  i^   product 
can  exhibit  a  randomly  fluctuating  phase  give  vanishing  contributions 
in  the  limit  n  ^  «>.  We  shall  now  identify  the  surviving  irreducible 
diagrams  in  the  ring  and  ladder  models  and  construct  the  corresponding 
closed  expressions  for  M,  (^  )• 

It  follows  from  (5.11)  and  (3-12)  that  the  first  three  of  relations 
(5.4)  are  satisfied  in  the  ring  model.   Thus  the  irreducible  diagrams 
Figs.  3b,  3c,  and  ^a  survive.   However,  the  1^   product  associated  with 
Fig.  4b  is 

^^         <^     ,       =  exp  i(9,   +9   )  +  i(0   +9.)  .  (5-7) 

"^kprs^^Brkp     ^j_  ^  ks    pr^    ^  sp    rk'J 

The  random  phase  of  this  product  does  not  vanish  when  p,  r,  and  s  are 
constrained  by  momentum  conservation  alone.   Consequently,  the  diagram 
does  not  survive.   The  further  irreducible  diagrams  which  do  survive 
in  the  ring  model  are  the  infinite  class  of  ring  diagrams,  whose  first 
two  members  are  shown  in  Fig.  7-   (if  the  incoming  and  outgoing  external 
lines  are  joined  together,  these  diagrams  form  symmetrical  rings.)   The 
/  products  associated  with  the  successive  ring  diagrams  have  the  phases 


(9,   +9   )  +  (9   +  9  ,  ,)  +  (9  ,  ,  +  9  .  )  . 
r   ks    pr    ^  rp    p'r'^    ^  r'p'    sk'J' 


1(9,   +9   )  +  (9   +  9  ,  ,)  +  (9  ,  ,  +  9  „  „)  +  (9^,,,  „  +  9   )  , 
[j  ks    pr^    ^  rp    p'r'    ^  r'p'    p  r  '    'r'p     sk  J 

(5.8) 


kl 


vhich  all  vanish,  by  (3.12). 

It  may  be  seen  from  the  rules  of  Sec.  k   that  the  ring  diagrams 
give  a  set  of  contributions  to  H  (^  )  which  resei±ile  a  geometric  series, 
They  may  be  siommed  easily  by  the  usual  methods  for  such  series.   If  we 
include  also  the  contributions  from  Figs.  5b,  5c,  and  ha,    the  final  re- 
sult for  Mj^(^j^)  is 


\(g  =  I^(\^\.^)\ 


+  P"'  II!  \-s\-k"(^c-a)Sp(^t)Ss(^c)\+p-s(^a+b-c^'       (^.9) 
psbc  -^  ^ 


where  V  '(C)  is  defined  by 

\'«d)  =  \to''   Ev,s^,(?,,)s^,.,(t,,.,)v,H?,).    (5.10) 


The  terms  in  (5-9)  proportional  to  N  arise  from  the  Hartree-Fock  diagrams 
Figs.  5b  and  5c.   To  verify  the  remainder  of  the  result,  we  may  solve 
(5-10)  by  iteration  to  yield  V  '(^^)  as  a  function  of  V  ,  and  then  substi- 
tute in  (5.9).   The  zeroth  term  V  in  the  series  for  V  ^{^^)   gives  rise 
to  the  contribution  of  Fig.  h-a.,    and  the  higher  terms  give  the  successive 
ring  diagram  contributions. 

Equation  (5.IO)  may  be  rewritten 


V  '(O  =  V 


l-^P"'EVqS   <^(?b'-d)Sp.(^b') 
p'b'  ^  -^   ^        ^ 


■1 

(5.11) 
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We  note  that  the  potential  components  V  for  different  values  of  q  are 
not  explicitly  mixed  in  this  expression^  a  fact  which  recalls  the  classi- 
cal ring  results  ohtained  in  Sec.  2.5'   We  may  call  V  '  ( ^ -, )  an  'effective 
potential'  for  the  ring  model.   In  field- theoretic  terminology,  it  is  a 
kind  of  vertex  operator.   From  a  particle  viewpoint,  we  may  consider 
k-s  and  ^    in  (5.9)  as  a  momentum  and  'energy'  which  are  taken  from  the 
incoming  particle,  transmitted  along  a  chain  of  intermediate  particles, 
and  finally  returned  to  the  original  particle  to  complete  the  ring. 

It  is  of  interest  to  compare  the  ring  summation  represented  by  (5*9) 
and  (5.10)  with  that  given  by  Montroll  and  Ward  .   The  latter  appears  to 
correspond  to  a  summation  over  ring  diagrams  for  the  primitive  linked- 
diagram  expansion.   It  omits  the  iterated  self -energy  corrections  which 
our  model  incJLudes. 

It  follows  from  (5-6)  and  (5-7)  that  all  of  the  relations  (5. if) 
are  satisfied  for  the  ladder  model.   Hence  the  irreducible  diagrams 
Figs.  5b,  5c,  ha,    and  ^b  all  survive.   Tlie  further  irreducible  diagrams 
which  survive  are  the  infinite  classes  of  ladder  diagrams  and  exchange 
ladder  diagrams.   The  first  two  members  of  each  class  are  shown  in  Figs. 
6  and  8,  respectively.   It  follows  from  (5-6),  (5.7)  and  the  rules  of 
Sec.  h   that  the  <j)   products  associated  with  corresponding  diagrams  in 
the  two  classes  have  the  same  value.   For  the  successive  diagrams  in 
either  class,  these  products  have  the  phases 


(-0  ,  +  e   )  +  (-9  +  e  ,  , )  +  (-0  ,  ,+  0  ,  )  , 
\2      pk    rs     ^   rs   r '  s '  '    ^   r '  s '    pk  J 


(-0  ,  +  9   )  +  (-9  +  9  ,  ,)  +  (-9  ,  ,+  0  „  „)  +  (-9  „  „+  9  ,  )L 
pk    rs       rs   r '  s '       r '  s '    r  s        r  s    pk  J 


(5.12) 
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which  vanish  identically.   The  ladder  and  exchange  ladder  diagrams 
exhaust  the  higher  irreducible  diagrams  which  survive  in  the  ladder 
model,   (it  should  be  noted  that  according  to  the  rules  of  Sec.  h   dia- 
grams like  Fig.  9a  and  Fig.  9b,  which  are  topologically  identical  with 
Figs.  6a  and  8b,  respectively,  are  not  to  be  counted  separately.) 
The  contributions  of  the  ladder  and  exchange  ladder  diagram 
sequences  can  easily  be  summed  in  closed  form,  in  a  similar  fashion 
to  the  ring  diagrams.   If  we  include  also  the  contributions  from  Figs. 
5b,  5c,  ^a,  and  4b,  the  final  result  for  the  ladder  model  is 


VU=  S  (^o^Vp^^ 


-^P"'  I]  (Vs^Vs^^srpk'^^a+b^V^b) 
psbc 


Ss(^c)Vp-s(^a+b-c)^  (5.13) 


where  V    '(f  ,  )  is  defined  by 
srpk  ^ ^a+b 


V_v'(^...J  =  V.-v  -  P'"  m  \-s'\'r-Dk'(^ 


1 
rpk  '^^a+b''  ~  "s-k   ^  ^rf—l,  's-s'"s'r'pk  ^^a+b' 


r  =  k  +  p  -  s  and 
with.r'  -  k  +  p  -  s'.   As  in  the  ring  model  result,  the  terms  in  (5-15) 

proportional  to  N  arise  from  Figs.  5b  and  5c.   If  (5-1^)  is  expanded 

by  iteration  and  the  result  for  V^  k'^^a+b^  ^^  substituted  in  (5-15), 

the  explicit  contributions  from  Figs,  ha.,    kh,    and  the  two  ladder  sequences 


kk  - 


are  obtained.   The  non-exchange  diagrams  all  arise  frcm  the  factor  V 

which  appears  in  (5-15)  and  the  exchange  diagrams  all  arise  from  the 

factor  +  V    when  this  procedure  is  carried  out.   The  quantities 
p-s 

V     ' ( C  ,  )  may  be  regarded  as  defining  an  effective  potential  for  the 
S2rpk  ^  ^a+b 

ladder  model. 

If  the  symmetry  constraint  (5-7)  is  relaxed  in  assigning  random 
values  to  the  phases  9  ,    an  abbreviated ladder  model  results  in  which 

K-P 

none  of  the  exchange  diagrams  survive.   The  abbreviated  model  violates 
the  second  of  conditions  (3.5),  and  in  consequence  it  requires  an  elabora- 
tion of  our  diagram  formalism:   The  two  junctions  which  make  up  a  vertex 
must  be  distinguished.   We  shall  not  discuss  this  model  further  here. 

Equations  (^+.12),  (5-9),  and  (5.10)  constitute  a  closed  set  of 
integral  equations  which  determine  S  (^  )  for  the  ring  model.   Similarly, 
(4.12),  (5.13),  and  (5.1^)  constitute  a  closed  set  for  the  ladder  model. 
These  equations  incorporate  extensive  classes  of  terms  from  the  primitive 
diagram  expansion  for  \i^g)    in  the  true  problem.   They  include  all 
primitive  diagrams  which  can  be  obtained  from  the  surviving  irreducible 
diagrams  by  repeatedly  inserting  these  irreducible  diagrams  into  them- 
selves and  each  other  as  self-energy  parts.   Examples  of  complicated 
primitive  diagrams  included  in  the  ring  and  ladder  models  are  shown  in 
Figs.  10a  and  10b,  respectively. 

6.    RELATION  BETWEEN  DISTINGUISHABLE  AND  INDISTINGUSIHABLE  PARTICLE 

MODELS 

We  wish  now  to  establish  a  correspondence  between  the  models  formu- 
lated in  Sees.  2  and  3  and  thereby  verify  that  the  thermodynamic  relations 
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obtained  in  Sec.  2  are  at  least  formal  classical  limits  for  the 
fermion  and  boson  models .   The  correspondence  is  already  suggested 
by  the  identity  of  the  lower  bounds  on  the  ring  and  ladder  model 
Hamiltonians  which  we  found  in  the  two  cases.   Our  procedure  here 
will  be  to  formulate  the  'iistinguishable-particle  problem  in  terms 
of  second-quantized  fields,  one  for  each  particle.   Then  we  shall 
appeal  to  two  assumptions:   The  thermodynamic  equivalence  of  canonical 
and  grand  canonical  ensembles  for  infinite  systems,  and  the  equivalence 
of  distinguishable  and  indistinguishable  particles  in  the  classical 
limit. 

The  Hamiltonian  of  a  system  of  N  distinguishable  but  similar 
particles  interacting  "Chrough  the  pair  potential  V(x)  may  be  written 

H=  Z]  Z]  Wn)\(n)  -^^i'  (^-^^ 

nm    kprs 

(n,m  =  1,2,  ...  ,  N). 

Here  we  have  introduced  a  separate  second- quantized  field  [labeled  (n)J 
for  each  particle,  and  we  restrict  the  system  to  those  states  which  are 
eigenstates  with  eigenvalue  one  for  all  the  nxomber  operators 

^     \(n)      =  E  \(n)'\(n)-  ^^'^^ 


^n)  =  4-"  \(n)   =  V«  ^k(n)  \{n) 
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The  commutation  relations  are 


LVn)'  S(m)  X  ^  ^™V'  ^^•''' 


The  restriction  of  each  field  to  one-particle  states  makes  the  choice 
of  plus  or  minus  commutator  immaterial. 

Let  us  consider  the  limit  N  — >  0°  with  p  and  p  =  n/q      constant.   We 
shall  assume  in  this  limit  that  the  canonical  ensemble  for  our  system 
gives  the  same  thermodynamics  as  a  grand  canonical  ensemble  in  which 
all  states  of  the  second- quantized  fields  are  admitted.   The  latter 
ensemble  is  chosen  so  that 


<^n)>  =^'  ^""^  =  5   <^"))  ^"^  '  ^^'^^ 

where 

<»(„,>  =-(=-^"-"'""(„,]/-[=-^"'-^'"}-  (fi-*) 

We  should  note  that  (6.5)  implies  |j  — ^  -  »  when  N  — $►  «  with  constant 
P  and  p.,  as  may  readily  be  verified  for  free  particles.   The  consequence 
is  that  the  one-particle  distribution  function  takes  the  Maxwell-Boltz- 
mann  form  in  the  limit,  as  it  must  for  consistency.   We  may  also  note 
that  the  variance  <(' (N/  ,^-1)  y   vanishes  in  the  limit. 
Now  let  us  replace  (6.2)  with  a  model  H.  , 

''^^^W-'S   VsVsVp,rH-s^k(n)\(m)\(m)'ls(n)^   ^^"^^ 


where  the  </'y.  are  the  same  parameters  as  in  Sec.  3.   The  thermodynamics 
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of  the  model  system  may  be  obtained  by  a  propagator  formalism  very 
similar  to  that  of  Sec.  h- .      Let  us  write 


8) 


and  then  define  S,  ( t  )  =  Y^  S,  ,  ^(t  )  m  terms  of  this  S,  (u,u')  by 
k^^a     '^—  k(n)^^a  k 

n    ^  ' 

(4.6).   Again,  it  will  not  affect  the  final  results  whether  the  fermion 
or  boson  case  is  taken.   With  these  definitions,  we  find  that  (h-.Q)    and 
(^.10)  hold  also  for  the  present  system,  with  N   -  ^H  ^.f  ) •   (O^  course, 
for  given  P  and  p,  the  chemical  potential  n  will  be  very  different  in  the 
present  case  than  in  Sec.  h,    as  we  have  noted  above.) 
Let  us  write 

Then  the  primitive  linked-diagram  expansion  for  S,  (^  )  is  given  by  the 
rules  of  Sec.  h   provided  the  following  changes  are  made: 

a)  Give  each  line  a  particle  label  as  well  as  a  momentum  label. 
Associate  with  each  line  a  factor  of  the  form  (6.9)- 

b)  With  each  vertex,  bearing  momentum  labels  p,  q,  r,  s  and  particle 
labels  n,  m,  m',  n'  as  shown  in  Fig.  11,  associate  a  factor 

^    pqrs  p-s  p+q,r+s  b+c,d+e  nn '  mm' 

c)  Sum  the  final  result  over  all  values  of  all  the  particle  labels, 
including  those  on  the  external  lines . 
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The  principal  difference  between  the  results  of  this  expansion  and 
that  of  Sec.  h   is  that  now  the  contributions  of  all  diagrams  with  particle 
exchange  vanish  in  the  liiTiit  N  — ^  °o,  Q   ^>  ^,    as  we  would  expect.   The 
formal  reason  this  happens  is  that  the  factors  of  the  form  6   , ,    severely 
restrict  the  summations  over  particle  labels  in  the  exchange  diagrams. 
In  the  diagrams  without  exchange,  each  closed  loop  represents  a  separate 
intermediate  particle  which  interacts  with  the  incoming  particle  or  with 
another  intermediate  particle.   Examples  are  shown  in  Fig.  12. 

Suppose  that  we  now  determine  the  ^  by  (5 -6)  and  (5-7)^  the 

conditions  for  the  fermlon  or' boson  ladder  model.   It  is  clear  from 
Sec.  5.3  that  if  similar  analysis  is  carried  out  for  the  present  case, 
the  surviving  diagrams  will  be  all  those  which  arise  from  the  irreducible 
diagrams  Figs.  5t>,  ^a,  and  the  ladder  sequence  illustrated  in  Fig.  6.   As 
we  have  just  noted,  there  are  no  exchange  contributions  for  the  present 
system,  and  consequently  the  exchange  ladder  sequence  of  Fig.  8  gives 
nothing. 

Let  us  compare  this  result  with  what  we  get  from  the  distinguishable- 
particle  ladder  model  of  Sec.  2.   We  replace  (6.7)  by  the  interaction 
Hamiltonian 

Hi  =  i  Hi  ^   \-s^n,m;k-s\+p,r+s\(n)\(m)\(m)%(n)^    ^^-^^^ 
nm   kprs        '    ' 


where  the  jzi        are  given  by  (2.9)  et  seq.   It  is  clear  that  the  only 
change  in  the  expansion  for  S,  ( ^  )  is  that  the  factor  '^^^Y^   ^    in  rule  b) 

K   a  pqx  o  P"-^ 

22 
above  must  be  replaced  by  the  factor  ^  V    .    We  readily  find  that 

■^       '^  n.m;p-s  p-s 
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precisely  the  same  diagrams  survive  in  the  present  case  as  for  the  pre- 
vious ladder  model.  We  shall  illustrate  the  equivalence  by  two  examples. 

Consider  first  the  diagram  of  Fig.  12c  and  equip  it  vith  momentum 
labels  as  in  Fig.  6a.   Its  contribution  contains  the  factors 


j-i[(k-s)+(s-s')+(s'-k)].d^^^l 


(6.11) 


according  to  our  rules  and  to  (2.9).   The  phase  of  this  expression  is 
identically  zero,  and  consequently  the  contribution  survives  when  it 
is  summed  over  the  particle  labels  n  and  m  in  accord  with  rule  c) . 
Next,  however,  consider  the  ring  diagram  Fig.  12d,  with  momentum  labels 
as  in  Fig.  7a- •   Its  contribution  contains  the  factors 

VVV  exp  -iq-(d    +d    +d   )L  (6.12) 

qqq     \_         ^    n,i  t ,m         m,nj 

where  q  =  k-s  and  we  have  used  the  momentum  conservation  relations.   The 
phase  of  this  expression  fluctuates  at  random  as  we  sum  over  all  values 
of  n,  t,    and  m,  so  that  the  contribution  does  not  survive  in  the  limit 
jj  ^  M^  il  -^  00  .   (The  contribution  from  the  special  value  q  =  0  also 
vanishes  in  the  limit.) 

The  equivalence  of  the  distinguishable  and  indistinguishable  ring, 
random- coupling,  and  Hartree-Fock  models  may  be  verified  in  a  similar 
fashion.  We  shall  give  one  further  illustration:  Consider  again  Fig.  12d, 
with  momentum  labels  as  in  Fig.  7a,  but  now  fix  the  i^^  ^^^   by  relation 
(2.19)  for  the  ring  model.   The  contribution  from  this  diagram  now  contains 
the  factors 
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V  V  V  exp 

q  q  q 


<i  (e   +0,   )  +  (e,   +9   )  +  (9   +9   m  , 


(6.1: 


and  we  see  from  (2.20)  that  the  phase  vanishes. 

We  have  seen  that  the  thermodynamics  of  our  second-quantized  dis- 
tinguishable-particle system  is  formally  the  same  in  the  limit  N  -^  «, 
0  — >  CO  whether  the  models  of  Sec.  2  or  those  of  Sec.  3  are  used.   Our 
argument  was  based  on  the  equivalence  of  canonical  and  grand  canonical 
ensembles  for  the  system,  and  also  implicitly  made  use  of  assumption  5 
of  Sec.  5-lj  which  underlies  all  our  work.   Now  let  us  make  the  further 
assimiption  that  in  the  classical  limit  (sufficiently  high  temperature 
and  low  density)  the  thermodynamics  of  fermion,  boson,  and  distinguish- 
able particle  systems  become  identical  in  the  true  problem  (all  jZ^'s  =  1), 
provided  the  same  values  of  P  and  (3  are  taken  in  each  case.   We  have  seen 
that,  except  for  exchange  diagrams,  the  models  of  Sec.  3  select  precisely 
the  same  diagrams  from  the  true-problem  expansion  for  S  (^  )  in  all 
three  cases:  fermion,  boson,  and  distinguishable-particle.   However,  the 
exchange  diagrams  do  not  contribute  in  any  event  in  the  classical  limit. 
It  then  follows  from  all  this  that  the  models  of  Sec.  2,  applied  to  a 
distinguishable  particle  system,  should  give  the  same  thermodynamics  in 
the  classical  limit  as  the  models  of  Sec.  3;  applied  to  fermlons  and 
bosons.   It  follows  that  the  classical  results  of  Sec.  2  for  A  should 
represent  classical  limits  of  the  ring,  ladder,  random- coupling,  and 
Hartree-Fock  models  for  fermions  or  bosons. 
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It  may  also  be  possible  to  investigate  the  relation  between  our 

two  kinds  of  models  by  using  the  formalisms  of  Montroll  and  Ward  or 

23 
Lee  and  Yang  ,  neither  of  which  require  second-quantization.   In  the 

method  of  Lee  and  Yang,  the  thermodynamics  of  dis-tinguishable  (Maxwell- 

Boltzmann)  particles  is  expressed  in  terms  of  a  'binary  kernel'  which 

is  determined  from  the  ordinary  two-particle  matrix  elements 

/k,piH.  |s,r"^  =  V,   5,        .  (6.lij-) 

\  •'■^'  i'  '  /     k-s  k+p,r+s  ^ 

Then  an  algorithm  is  used  to  obtain  results  for  fermions  or  bosons. 
To  attempt  an  expression  of  our  models  in  this  formalism,  we  woiild 
replace  (6.1^)  by  the  model  matrix  element 

<Ck,p|H.  |s,r\     =  i2$    ,   V,   8,  (6.15) 

\  '-^1  ii  '   /  n,m   '^n,m;k-s  k-s  k+p,r+s  *■ 

for  the  Sec.  2  models,  or  by 

<'k,p|H.  Is.rN     =  y<,    V,   5,  (6.l6) 

\  '-^1  i'  -  /  n,m   "^kprs  k-s  k+p,r+s  ^ 

for  the  Sec.  5  models.   Here  n  amd  m  denote  the  pair  of  particles  for 
which  the  matrix  element  is  evaluated.  We  have  not  explored  this  pro- 
cedure . 

Recognition  of  the  thermodynamic  results  of  Sec.  2  as  classical 
limits  for  the  fermion  and  boson  models  may  provide  some  useful  insights 
into  the  behavior  of  the  latter.   For  example,  the  lack  of  saturation 
found  in  the  classical  ladder  model  suggests  that  a  similar  lack  charac- 
terizes the  fermion  and  boson  ladder  models.  However,  the  correspondence 
between  the  Sec.  2  and  Sec.  3  models  also  leads  to  a  rather  discouraging 
general  observation.   It  points  out  how  modest  are  the  presently  feasible 
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quantum-mechanical  diagram  summations  compared  to  known  classical  ones. 
Our  fermion  and  boson  ladder  and  ring  summations,  carried  out  in  Sec. 
5.5,  are  more  comprehensive  than  those  usually  employed;  they  include 
infinite  classes  of  self-energy  corrections  which  usually  are  omitted. 
Nevertheless,  we  have  seen  that  the  ladder  summation  corresponds  in  the 
classical  limit  to  just  one  term  in  the  Mayer  irreducible  cluster  expan- 
sion.  Similarly,  our  ring  summation  represents  only  a  partical  contri- 
bution from  each  term  of  the  classical  irreducible  ring  diagram  sequence. 
The  familiar  Montroll-Mayer  summetion  (2.29)  appears  to  correspond,  in 
the  quantum-mechanical  case,  to  retaining  both  ladder  and  ring  summations, 
together  with  the  self-energy  corrections  of  all  orders  obtained  by 
repeatedly  inserting  the  retained  irreducible  diagrams  as  self-energy 
parts . 

7.    DISCUSSION 

In  the  present  paper  we  have  obtained  formally  exact  closed  equations 
which  express  the  statistical  mechanics  of  a  class  of  model  Hamiltonians 
for  infinite  many-body  systems.   The  potential  value  of  these  equations 
lies  largely  in  the  fact  that  certain  of  the  models,  the  ring  and  ladder 
models,  share  important  boundedness  properties  with  the  true  many -body 
Hamiltonian:  The  eigenvalues  of  the  ladder-model  Hamiltonian  are  non- 
negative  if  the  pair  potential  V(x)  is  purely  repulsive,  and  those  of  the 
ring-model  Hamiltonian  are  bounded  from  below  if  V(x)  is  bounded  and  has 
a  non-negative  Fourier  transform. 
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These  properties  suggest  that  the  ring  and  ladder  models,  with 
appropriate  V(x),  should  have  a  meaningful  statistical  mechanics  even 
in  the  zero-temperature  limit.  Moreover,  the  structure  of  the  models 
is  such  that  they  embody  some  important  qualitative  dynamical  features 
of  the  time  many-body  system.   For  example,  we  may  expect  that  dissi- 
patlve  damping  of  single-particle  excitations  survives  in  the  models. 
This  is  particularly  clear  in  the  distinguishable-particle  formulation 
of  the  models:  each  particle  interacts  individually  with  every  other 
particle.  A  similar  situation  exists  for  the  fermion  and  boson  models, 
but  there  it  is  more  natural  to  think  of  interaction  among  individual 
momentiim  modes  rather  than  among  individual  particles. 

The  remarks  just  made  suggest  that  the  ladder  or  ring  models  may 
be  appropriate  for  investigating  whether  the  sharp  Fermi  surface  of  an 
infinite  system  of  uncoupled  fermions  at  zero  temperature  persists 
when  the  particles  are  coupled.   Similarly,  the  model  solutions  may  be 
of  aid  in  deciding  whether  singular  occupancy  of  the  zero-momentum  state, 
which  characterizes  an  infinite  free-boson  system  at  very  low  temperatures, 
persists  when  the  particles  are  coupled.   (We  anticipate  here  the  extension 
of  our  analysis  to  low- temperature  boson  systems  which  will  be  carried  out 
in  the  following  paper.)   The  natural  way  in  which  the  effective  density 
components  p  appear  in  the  ring  model  suggests  that  it  may  be  appropriate 
for  investigating  phonon-like  excitations  and  other  collective  phenomena. 

However,  any  confidence  which  our  rigorously  bounded  model  Hamil- 
tonians  may  inspire  in  a  given  problem  does  not  automatically  extend  to 
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the  closed  equations  which  we  have  derived  for  the  model  propagators. 
First  we  must  establish  that  these  equations  are  exact  descriptions  of 
the  models  in  actuality,  as  well  as  formally.   This  has  not  been  done 
in  the  present  paper.   As  we  discussed  in  Sec.  ").!,    a  fundamental  con- 
vergence question  relating  to  extremely  high-order  diagrams  is  involved. 
We  shall  state  this  question  more  precisely  in  Sec.  h   of  the  following 
paper,  using  generalized  stochastic  models  which  yield  our  formally 
closed  equations  for  finite  as  well  as  infinite  systems.   In  Sec.  7  of 
that  paper,  we  shall  outline  what  we  hope  is  the  basis  for  a  satisfying 
justification  of  our  closed  model  propagator  equations. 

The  analysis  in  the  following  paper  is  concerned  almost  exclusively 
with  the  indistinguishable-particle  models,  and  we  shall  not  attempt 
there  to  offer  explicit  justification  for  the  assumptions  which  underlie 
the  classical  results  of  the  present  Section  2.   However,  we  have  already 
found  a  degree  of  support  for  these  results:   The  closed  expressions  for 
the  Helmholtz  free  energy  obtained  in  Sec.  2  reproduced  precisely  the 
rigorous  lower  bounds  on  the  ring  and  ladder  model  potential  energies . 

Our  third  model,  the  random- coupling  model,  exhibited  no  lower  bound 
on  the  potential  energy  per  particle  in  the  limit  of  an  infinite  system. 
The  results  of  Sec.  2.k   suggest  that  this  is  the  case  whatever  the  form 
of  V(x) .   Consequently,  we  must  expect  the  random- coupling  model  not  to 

give  sensible  statistical-mechanical  results  at  zero  temperature  in  either 

2k 
the  classical  or  the  quantum-mechanical  case. 

The  assured  failure  of  the  random- coupling  model  at  zero  temperature 

may  point  a  moral.   We  have  seen  in  Sec.  5  that  this  model  corresponds 

to  taking  just  the  lowest  few  diagrams  in  the  irreducible  diagram  expansion 
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of  the  propagator  for  the  true  many-body  Hamiltonian.   The  failure  of 
the  model  suggests  that  a  term-by-term  treatment  of  the  irreducible 
diagram  expansion  for  the  true  problem  may  not  be  justified  at  very  low 
temperatures,  ajid  that  results  obtained  from  such  a  treatment  shoiild  be 
viewed  with  caution. 

A  general  question  raised  by  the  present  paper  is  whether  there 
exists  an  infinite  sequence  of  stochastic  models  which  correspond  to 
more  and  more  comprehensive  (but  siimmable)  classes  of  terms  from  the 
linked  diagram  expansion  for  the  true  many-body  problem.   We  have  so 
far  not  succeeded  in  constructing  substantially  more  elaborate  models 
than  those  presented  here.  An  obvious  next  step  is  to  seek  a  model  that 
combines  both  ring  and  ladder  summations  so  as  to  correspond  to  the 
classical  Montroll-Mayer  ring-diagram  summation.   On  the  basis  of  a 
preliminary  investigation,  we  offer  the  following  opinion:   If  such  a 
model  can  be  constructed  within  the  general  formal  framework  of  Sec. 
3.1,  it  probably  can  be  achieved  only  by  allowing  the  parameters  i^j^  ^ 
to  have  stochastically  distributed  moduli  as  well  as  phases. 
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APPEITOIX 

We  wish  to  discuss  here  two  topics  which  arose  in  Sec.  2:  the 
derivation  of  (2.11)  for  the  classical  models^  and  the  significance 
of  the  condition  (Bp/Sp)   <  0  for  the  models 

Equation  (2.11)  for  the  true  problem  usually  is  derived  on  the 
basis  of  the  grand  canonical  ensemble  and  rigid-wall  boundary  condi- 
tions.  It  then  is  formally  exact  in  the  limit  0  — >».   [see,  for 
example^  T.L.  Hill,  Statistical  Mechanics  (McGraw-Hill  Book  Company, 
Inc.,  New  York,  1956),  chap.  5.]   However,  with  our  cyclic  V(x)  the 
factoring  of  reducible  cluster  integrals  into  irreducible  integrals 
is  exact  for  any  fi.   Consequently,  (2.11)  is  formally  valid  for  any 
n,  with  the  grand  canonical  ensemble.   In  the  limit  N  — >  o^,  Q   — ^  oo 
it  is  usually  considered  immaterial  whether  the  grand  canonical  or 
canonical  distribution  is  used,  and  we  shall  assiome  that  this  is  so 
here. 

In  the  case  of  our  models,  a  grand  canonical  ensemble  can  be 
formed  by  considering  each  of  the  N  particles  in  the  canonical 
ensemble  as  a  separate  species  and  taking  an  accivity  such  that  the 
mean  total  number  of  particles  over  the  grand  ensemble  is  N.   [in 
doing  this,  we  may  extend  (2.5)  -  (2.8)  to  include  the  case  n  =  m, 
so  as  to  allow  interaction  among  particles  of  the  same  species.]   The 
derivation  of  (2.11)  in  the  limit  N  — >  00^  i]  —^^,    with  B  interpreted 
as  in  the  text,  then  depends  upon  two  facts .   First,  a  reducible  cluster 
integral  involving  any  given  a  +  1  particles  factors  exactly  into  irre- 
ducible cluster  integrals.   As  in  the  true  problem,  this  is  true  for 
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any  Q.  because  V^''"(x)  is  cylcic.   Second,  if  B  is  averaged  over  all  N 

choices  of  the  species  of  any  one  of  the  particles  in  the  cliister,  then 

(in  the  limit  N  -^  «.  ,  H  ->  oo  )  the  result  is  the  same  for  all  choices 

of  the  species  of  the  other  particles,  except  for  a  set  of  choices  of 

relative  measure  zero.   This  can  be  seen  for  each  of  the  models  by 

analysis  similar  to  that  used  in  the  text  to  evaluate  the  B  .  We  shall 

a 

not  give  here  a  detailed  derivation  of  (2.11)  for  the  models.   It  is 
straightforward  once  the  two  facts  just  stated  are  established. 

As  in  the  true  problem,  we  assume  in  the  text  that  the  final  averages 
for  an  infinite  system  are  independent  of  whether  the  grand  canonical  or 
canonical  ensemble  is  used.   Tlie  canonical  ensemble  is  taken  in  Sec.  2 
because  it  makes  the  discussion  simpler.   In  the  quantum-mechanical  treat- 
ment of  Sec.  6,  the  grand  ensemble  is  employed. 

The  derivation  and  analysis  of  (2.11)  for  the  models  takes  a  much 
more  elegant  form  if  one  uses  the  generalized  models  described  in 
Appendix  A  of  reference  11.   Then,  with  the  grand  ensemble,  the  formula- 
tion in  terms  of  averaged  irreducible  clusters  is  exact  for  all  N  and  0. 

The  statement  in  the  text  that  (Sp/Sp)   <  0  implies  instability 
must  be  carefully  qualified.   Actually,  the  condition  (Sp/Sp)   <  0,  or 
even  p  <  0,  does  not  necessarily  imply  instability  for  any  of  our  models, 
if  p  is  defined  as  p  (Sa/Sp)  .   This  is  because  of  the  peculiar  way  in 
which  the  potentials  V^'"'(x)  are  constructed.   The  simplest  illustration 
is  provided  by  the  Hartree-Fock  model  as  described  by  (2.^)   and  (2.35). 
Suppose  that  we  have  V  <  0,  so  that  the  potential  energy  per  particle 
(which  is  just  A-A  for  this  model)  decreases  without  limit  as  p  increases. 
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that  is^  as  we  pack  more  particles  into  a  fixed  volume  H.   By  (2.55)? 
we  have  p  <  0  and  (3p/Sp)   <  0,  if  P  is  high  enough.   However,  since 
each  particle  moves  in  a  uniform  potential,  the  potential  energy  is 
independent  of  configuration,  for  given  R  and  ^,    and  the  system  is  not 
unstable.   Instability  arises  only  if,  as  we  pack  the  N   particles  into 
a  closer  configuration,  we  decrease  correspondingly  the  volume  of  the 

^     ,    the  potential  energy  becomes  in- 


creasingly negative  and  the  system  can  collapse  catastropliically. 

It  is  likely  that  the  relation  (Bp/Bp)   <  0  has  a  similar  inter- 
pretation for  all  the  models  treated  in  this  paper.   This  condition  need 
not  imply  instability  if  H  is  fixed,  which  it  must  be  as  we  have  defined 
the  models.   However,  if  we  were  to  allow  fi  to  vary  in  accordance  with 
the  actual  gross  volume  occupied  by  the  particles,  there  would  be  in- 
stability.  Having  H  vary  would  actually  be  a  physically  appropriate 
procedure,  as  the  Hartree-Fock  example  suggests.   For  this  reason,  we 
consider  (Sp/9p)   <  0  to  be  an  instability  indication  in  making  a 
physical  interpretation  of  our  models. 

It  should  be  noted  also  that  the  lack  of  a  lower  bound  on  the 
potential  energy  of  a  system  does  not  preclude  a  stable  thermodynamics. 
If  the  density-of-states  G'(e)   (where  E  is  the  total  energy)   decreases 
faster  than  exponentially  as  E  — >  -  oo  ,  a  stable  thermal  equilibrium 
can  exist  for  all  finite  (3.   If  the  decrease  is  slower  than  exponential, 
equilibrium  will  be  impossible  at  any  p,  while,  if  the  decrease  is  expo- 
nential, eqiillibrium  can  exist  only  if  p  is  less  than  a  critical  value. 
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We  wish  finally  to  note  the  conclusion  of  L.  van  Hove  [physica  I5, 
951  (19^9)]  that  (Sp/Sp)   <  0  cannot  be  an  exact  theoretical  result  for 
a  gas  of  pair-wise  interacting  particles.  Althovigh  it  seems  assuredly 
valid  for  actual  physical  systems,  van  Hove's  result  does  not  appear  to 
be  applicable  here.   Our  models  violate  an  assumption  basic  to  his 
analysis:   The  v  '  (x)  are  defined  in  such  a  way  that  it  is  not  possible 
to  divide  the  system  into  effectively  non-interacting  macroscopic  sub- 
volumes.   (Van  Hove's  result  clearly  is  not  valid  for  the  Hartree-Fock 
example  discussed  above.)   It  should  be  stressed,  however,  that  the 
results  obtained  in  Sec.  2  are  not  rigorously  justified  by  the  analysis 
we  have  presented.   As  we  have  noted,  there  is  a  fiindamental  convergence 
question  involved.   We  shall  take  up  this  question  in  the  following 
paper. 
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1.  One  way  to  obtain  this  result  is  the  following.   Take  X  <  0  (pure 
attractive  potential) .   Then  every  diagram  in  the  irreducible  cluster 
expansion  for  pressure  gives  a  negative  contribution.   The  total  number 
of  diagrams  of  order  n  increases  with  n  faster  than  any  exponential^  and 

a  consequence  is  that  the  pressure  comes  out  negatively  infinite  no  matter 
how  small  |x|  is.   (Physically,  this  means  the  system  will  collapse.)   On 
the  other  hand,  for  X  >  0,  the  pressure  must  approach  the  perfect  gas 
value  as  |x|  ^0.   Therefore  the  pressure  is  a  nonanalytic  function  of 
X  at  X  =  0.   [Cf.  T.D.  Lee  and  CN.  Yang,  Phys .  Rev.  105,  1119  (1957)-] 
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6.  A.  Abrikosov,  L.P.  Gor'kov,  and  I.E.  Dzyaloshinskii,  J.  Exptl.  Theoret. 
Phys.  (U.S.S.R.)  56,  900  (1959)  [translation:  Soviet  Phys.  JETP  9^  636 
(1959)]. 

7.  P.O.  Martin  and  J.  Schwinger,  Phys.  Rev.  II5,  15^2  (1959). 

8.  J.M.  Luttinger  and  J.C.  Ward,  Phys.  Rev.  II8,  1^17  (I960). 

9.  D.N.  Zubarev,  Usp.  Fiz.  Nauk.  71,  71  (1960)  [translation:  Soviet  Phys. 
Uspekhi  3,    320  (1960)] . 

10.   Extensive  further  bibliography  is  given  in  references  7  and  9. 
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11.  R.H.  Kraichnan,  Report  HT-10,  Division  of  Electromagnetic  Research, 
Institute  of  Mathematical  Sciences,  New  York  University  (I96I,  to  appear). 

12.  See  the  Appendix. 

15'   As  (3  —>  »  ,  we  have  A  -»  0,  and  at  zero  temperature  A  becomes  just 

the  potential  energy  per  particle. 

1^.   The  following  may  make  clear  how  this  can  happen.   Take  n   finite 

(but  »  r  )  and  place  an  arbitrarily  large  number  of  particles  into  the 

cube  in  any  desired  positions  x  .  Whatever  the  number  of  particles,  and 

whatever  their  positions,  it  is  clear  that  for  every  pair  n,  m  there  will 

be  many  possible  choices  of  d    such  that  Ix  -  x  -  d   I  >  r  . 
n,m  'nmn,m'0 

15.  This  statement  must  be  carefully  qualified.   See  the  Appendix  for 
a  discussion  of  the  condition  (Sp/3p)   <  0. 

16.  D.  Pines  and  D.  Bohm,  Phys.  Rev.  85,  358  (1952). 

17-   E.W.  Montroll  and  J.E.  Mayer,  J.  Chem.  Phys.  9,  626  (19^1). 

18.  Our  result  also  resembles  (except  for  the  terra  p  PW  )  the  classical 
limit  of  a  quantum- me chani cal  ring  summation  given  in  reference  5  [Eq. 
(5-1^)  of  that  ref erencej .   However,  that  summation  is  based  on  an  activity 
rather  than  a  density  expansion.   It  corresponds  to  a  sum  of  pure  ring  dia- 
grams from  the  primitive  (reducible)  classical  cluster  expansion. 
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20.  See  N.M.  Hugenholtz  and  D.  Pines,  Phys.  Rev.  II6,  ^4-89  (1959)- 
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21.   If  the  vertices  and  jtmctions  were  considered  distinguishable,  one 
would  obtain,  for  each  of  our  distinct  diagrams,  a  total  of  2  nl  diagrams, 
all  of  which  would  give  identical  contributions  to  S  (^  ).   If  one  counts 
these  diagrams  separately,  which  we  do  not  do,  the  contribution  per  diagram 
obtained  by  o\ir  irules  1-7  must  be  multiplied  by  l/2  nl   The  latter  procedure 
is  adopted  in  reference  8. 


22.  We  may  formally  extend  the  definition  of  the  p   -v  ^°  include  the 
case  n  =  m,  as  we  did  in  Sec.  2.5.   This  case  represents  a  vanishing  con- 
tribution here,  in  the  limit  N  -^  «. 


23.  T.D.  Lee  and  CN.  Yang,  Phys .  Rev.  11_3;>  II65  (1959);  Phys .  Rev.  II7, 
22  (i960). 

24.  This  does  not  necessarily  exclude  the  random- coupling  model  for  con- 
densed boson  systems.   The  classical  equation  of  state  (2.53)  suggests 
admissibility  at  any  given  p  and  p  provided  V(x)  is  weak  enough.   Thus  the 
model  may  be  admissible  quantum-mechanically  even  below  the  \-point,  if 
V(x)  is  weak  enough.   We  regard  this  argument  with  strong  suspicion. 
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(a) 


(c) 


Fig.  1.   Some  irreducible  cluster  diagrams. 


(a) 


(b) 


Fig.  2.   (a)  A  vertex.   (b)  A  labeled  vertex. 
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(a: 


(b) 


(c) 


Fig.  5.   The  zeroth  and  first-order  diagrams  for  S  (0, 
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(Q) 


(b) 


Fig.  k.      The  second-order  irreducible  diagrams  for  S^(?g. 


65  - 


Fig.  5.   A  primitive  diagram  which  survives  in  the  Hartree-Fock 


(a) 


(b) 


Fig.  6.   The  third- order  and  fourth-order  ladder  diagrams. 
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(b) 


Fig.    7-      The   third-order  and  fourth-order  ring  diagrs 
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(a) 


(b) 


Fig.  8.   The  third-order  and  fourth-order  exchange  ladder  diagrams. 
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(a) 


(b: 


Fig.  9.      (a)  A  diagram  redundant  with  Fig-  6a, 
(b)  A  diagram  redundant  with  Fig.  8b. 
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(a) 


(b: 


Fig.  10.   Examples  of  higher  primitive  diagrams  contributing 
in  the  ring  model  (a)  and  the  laddex-  model  (b). 
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Fig.  11.   A  labeled  vertex  for  the  second-quantized 
distinguishable  particle  system. 
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Fig.  12.   Some  diagrams  for  the  system  of  distinguishable  particles. 
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